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ABSTRACT 

We analyzed absorption features arising from interstellar neutral carbon that appeared in the 
UV spectra of 89 stars recorded in the highest resolution echelle modes of the Space Telescope 
Imaging Spectrograph on HST so that we could determine the relative populations of collision- 
ally excited fine-structure levels in the atom's electronic ground state. From this information, 
in combination with molecular hydrogen rotation temperatures, we derive the distribution of 
thermal pressures in the diffuse, cold neutral medium. We find a lognormal pressure distribution 
(weighted by mass) with a mean in \og(p/k) equal to 3.58 and an rms dispersion of at least 
0.175 dex that plausibly arises from turbulence with a characteristic Mach number in the range 
1 < M < 4. The extreme tails in the distribution are above the lognormal function however. 
Overall, pressures are well correlated with local starlight intensities and extreme kinematics, and 
they show some anticorrelation with kinetic temperatures. A subsample restricted to low am- 
bient UV intensities reveals a mode in the distribution of log(p/fc) that is nearly the same as 
the complete sample, but with a strong negative skewness created by a near absence of a tail at 
high pressures. Approximately 23% of this gas is at a pressure that is below that allowed for 
a static cold neutral medium. Accompanying nearly all of the gas is a small fraction (~0.05%) 
that has an extraordinarily large pressure, log(p/fe) > 5.5, and this condition is more prevalent 
at high velocities or for regions with enhanced starlight densities. This survey suggests that the 
dispersion of thermal pressures in the cold, neutral ISM is predominantly governed by micro- 
scopic turbulence driven by star-forming regions, with some additional effects from macroscopic 
events (e.g., SN explosions), and these measurements provide constraints for future studies of the 
broader impact of turbulence on the ISM and star formation. 

Subject headings: ISM: atoms - ISM: kinematics and dynamics - ISM: lines and bands - techniques: 
spectroscopic - turbulence - ultraviolet: ISM 
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1. Introduction 

In his commentary on the classic paper by 
Chandrasekhar & Munch (1952) on brightness 
fluctuations in the Milky Way, Scalo (1999) pre- 
sented an insightful discussion of a dichotomy in 
our perception of the structure of the diffuse inter- 
stellar medium (ISM) of our Galaxy. On the one 
hand, we might view the ISM in terms of a collec- 
tion of isolated, dense clouds enveloped in a more 
tenuous medium, a concept that has directed our 
thinking on the establishment of discrete "phases" 
of the ISM with well established spatial domains 
and vastly different properties that can be justi- 
fied on some fundamental physical grounds (Field 
et al. 1969; McKee & Ostriker 1977; Burkert & 
Lin 2000; Vazquez-Semadeni et al. 2000; Bran- 
denburg et al. 2007). On the other hand, much 
of the ISM can be viewed as a continuous fluid 
medium containing a texture of seemingly random 
fluctuations in density, velocity and temperature. 
Adherents to this second picture [e.g., Ballesteros- 
Paredes et al. (1999)] view clouds as an illu- 
sion created by the most extreme fluctuations in a 
turbulent medium having an extraordinarily high 
Reynolds number. While this is undoubtedly true, 
we must also acknowledge the presence of nearly 
static, sharp boundaries between different media, 
as revealed by dark clouds with well defined edges 
that have been sculpted by ionization, dissocia- 
tion, and evaporation/condensation fronts. Both 
pictures have their utility in exploring important 
issues on the multitude of processes that can occur 
within the ISM. 

As we switch our perspective from morpholog- 
ical to dynamical properties of the ISM, we find 
that over macroscopic scales the motions of gases 
in our Galaxy can be governed by the injection 
and dissipation of mechanical energy from a wide 
range of energy sources that include supernova ex- 
plosions (McKee & Ostriker 1977; McCray & Snow 
1979; Mac Low et al. 1989; Kim et al. 2001; dc 
Avillez & Breitschwerdt 2005a), disturbances from 
newly formed H II regions (Lasker 1967; Tenorio- 
Tagle 1979; Rodriguez-Gaspar & Tenorio-Tagle 
1998; Peters et al. 2008), stellar mass loss (Abbott 
1982; McKee et al. 1984; Owocki 1999), infalling 
gas clouds from the Galactic halo (Wakker & van 
Woerden 1997; Santillan et al. 1999, 2007) , bipo- 
lar jets from star forming regions (Bally 2007), 



shocks in spiral arm density waves (Roberts et 
al. 1975), and the magnetorotational instabililty 
driven by differential galactic rotation (Pinotek & 
Ostriker 2004) . These processes play a strong role 
in creating recognizable, discrete structures and 
flows of material in the ISM, but ultimately some 
of the energy from the resulting compressions and 
vorticity will also be transformed into random tur- 
bulent motions. Transient structures of small sizes 
can arise from the cascade of larger turbulent cells 
into small ones or be created in the interface re- 
gions between colliding gas flows (Audit & Hen- 
nebelle 2005). Turbulence can also be fed by in- 
stabilities in phase transition layers (Inoue et al. 

2006) or the weak driving forces that arise from 
the thermal instability of the ISM (Kritsuk & Nor- 
man 2002b; Koyama & Inutsuka 2006), the latter 
of which can be sustained by abrupt changes in 
the heating rate from UV radiation (Kritsuk & 
Norman 2002a). 

Over the past several decades, much progress 
has been made in the study of magnetohydrodymi- 
cal (MHD) turbulence in the ISM. Our under- 
standing of this phenomenon has been facilitated 
by the rapid emergence of powerful 3-dimensional 
computer simulations, and its existence is sup- 
ported by observations of column density distribu- 
tions, velocity statistics, cloud morphologies, de- 
viations in magnetic fields, and various kinds of 
disturbances in the propagation of radio waves 
in ionized media [for a review, see Elmegreen & 
Scalo (2004)]. This phenomenon is influential on 
the heating, chemical mixing, radio wave propaga- 
tion, and cosmic ray scattering in the ISM (Scalo 
& Elmegreen 2004). Within denser environments, 
turbulent processes are expected to have a strong 
influence on the fragmentation of density concen- 
trations just before and during the earliest stages 
of gravitational collapse that leads to star forma- 
tion (Mac Low & Klessen 2004; McKee & Ostriker 

2007) . 

Both the coherent dynamical phenomena and 
the smaller scale turbulent motions have an in- 
fluence on pressures in the ISM. These pressures 
appear in many forms: thermal, magnetic, dy- 
namical, and the indirect effects of cosmic rays, 
and their collective magnitude amounts to about 
p/k = 2.5 x 10 4 cm" 3 K0 that is established by 



Throughout this paper, we quantify pressures in terms of 
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the hydrostatic equilibrium of gaseous material in 
the gravitational potential of the Galactic plane 
(Boulares & Cox 1990). Except for very hot me- 
dia (T > 10 5 K) that have been created by shock 
heating from supernova blast waves (Cox & Smith 
1974; de Avillez & Breitschwerdt 2005b) or that 
reside within wind-blown bubbles around stars 
(Castor et al. 1975; Weaver et al. 1977), the ther- 
mal pressures of the ISM represent a small fraction 
(about one-tenth for T ~ 100 K) of the total pres- 
sure (de Avillez & Breitschwerdt 2005a). While 
thermal pressures and their variability may seem 
unimportant in the dynamical development of the 
general ISM, they nevertheless can be influenced 
by stochastic dynamical effects and thus can pro- 
vide us with useful information. 

In this paper, we make use of the fact that the 
two upper fine-structure levels in the electronic 
ground state of the neutral carbon atom, with ex- 
citation energies E/k = 23.6 and 62.4 K, are eas- 
ily excited and de-excited by collisions with neu- 
tral and charged particles at typical densities and 
temperatures within the diffuse, cold gas in the 
Galactic plane. The balance between the effects of 
these collisions and spontaneous radiative decays 
(at wavelengths 609 and 371 /im) establishes fine- 
structure level population ratios that can serve as 
an indicator of the local density and temperature 
of the C I-bearing material. We sense these popu- 
lation ratios by observing the UV multiplets of C I 
that appear as foreground absorption features in 
the spectra of hot stars recorded at high spectral 
resolution. 

The first widespread study of C I fine-structure 
excitations was carried out by Jenkins & Shaya 
(1979), who analyzed observations that came from 
the UV spectrograph on the Copernicus satellite. 
Those observations and a more comprehensive sur- 
vey by Jenkins et al. (1983) were primitive by 
today's standards for observing C I features set 
by spectrographs on the Hubble Space Telescope 
(HST) (Smith et al. 1991; Jenkins et al. 1998; 
Jenkins 2002), but they nevertheless established 
an early framework for determining the average 
thermal pressures in the ISM and their variations 
from one location to the next. In addition to 



p/k in the units cm -3 K instead of simply p in the units 
of dyne cm" 2 or erg cm -3 . Our representation facilitates 
comparisons with actual densities and temperatures in the 
ISM. 



these general surveys, special studies sensed ex- 
treme positive deviations in pressure from the C I 
features in spectra of stars within and behind the 
Vela supernova remnant (Jenkins et al. 1981, 
1984, 1998; Jenkins & Wallerstein 1995; Waller- 
stein et al. 1995; Nichols & Slavin 2004), indicat- 
ing that the blast wave overtook and compressecH 
small clouds in the medium that surrounded the 
explosion site (Chevalier 1977). 

A new advance in the study of C I excitation in 
the general ISM arose from the study by Jenkins & 
Tripp (2001) (hereafter JT01), who used the high- 
est resolution configurations of the echelle spectro- 
graph in the Space Telescope Imaging Spectrograph 
(STIS) (Kimble et al. 1998; Woodgate et al. 1998) 
to observe C I in the spectra of 21 stars. An impor- 
tant breakthrough in this work was to make use 
of the ability of STIS to record many different C I 
multiplets simultaneously, which allowed JT01 to 
benefit from a special analysis technique that they 
developed to unravel the blended absorption pro- 
files arising from the three fine-structure levels. 

The current study of thermal pressures expands 
on the work of JT01, once again using the analysis 
method employed earlier, but with some technical 
improvements outlined in Appendix [AJ In Sj2] we 
describe our new coverage of sightlines that sig- 
nificantly expands on the limited selection of stars 
that were studied by JT01, but we caution that 
one must be aware of a few, mostly unavoidable, 
selection biases in the sampling. In <J3] we review 
the basic principles of the analysis, but leave it 
to the reader to consult JT01 for a more detailed 
description of the mathematical method. This sec- 
tion also introduces our fundamental approach to 
interpreting the population ratios in terms of the 
local density and temperature of the C I-bearing 
gas, a method originally developed by Jenkins & 
Shaya (1979). Table [2] in this section lists the 89 
target stars whose spectra were analyzed in this 
study, along with some relevant information about 
the foreground regions probed by the sightlines. 
For each sightline, we show in Table [3] some com- 
posite interstellar conditions that we derived and 
some reflections on their significance in $7] In Ta- 
ble [4] we provide more detailed, machine- readable 



2 We note that Nichols & Slavin (2004) proposed some pos- 
sible alternative explanations for producing an excess of 
excited C I. 



3 



information for each velocity bin in the entire sur- 
vey. Section |4] describes a number of additional 
factors that require some consideration when the 
C I results are interpreted, including ways to es- 
timate the total amount of gas that accompanies 
the C I (<2jJ, the kinetic temperature of the gas 
( §4.2p . the mix of particles that can excite the up- 
per levels ( ^4.3|) . and the local intensity of starlight 
( §4.4p . The starlight intensity must be known in 
order to correct for optical pumping of the fine- 
structure levels (as described in Sj6j, and, as we 
demonstrate in £19. 11 such intensities seem to be 
correlated with pressure. 

Section [5] discusses how we interpret our find- 
ing that virtually all of the measurements do not 
agree with the expected fine-structure populations 
for collisional excitation for any uniform values 
of local density or temperature. Here, we intro- 
duce the idea that small admixtures of gas at ex- 
traordinarily high pressures accompany virtually 
all of the gas at normal pressures, and in ^10.5.11 
we examine (and ultimately reject) some alterna- 
tive explanations of the observed deviations. Some 
discussions on the implications and possible ori- 
gins of the high pressure component appear in 
^10.5.2| - [TUX1 and we show in SjS] and t JlO.5.41 
that the gas fractions at high pressures are ac- 
centuated in material that is moving rapidly or is 
exposed to a high intensity of starlight. 

In Sj5]we derive the distribution function for the 
mass-weighted thermal pressures in the dominant 
low pressure regime for two samples: (1) all of 
the gas and (2) gas that is well removed from 
intense sources of UV radiation. For the conve- 
nience of those who wish to compare our results 
with computer simulations of turbulence based 
on volume-weighted distributions, we convert our 
mass- weighted sampling to a volume- weighted one 
in §9.2[ but with the precarious assumption that 
the gas responds to pressure perturbations with a 
single polytropic index 7 (which is left as a free pa- 
rameter). In we address the possibility that 
the width of the pressure distribution understates 
the true dispersion of pressures in the ISM, due to 
the fact that we view in each velocity bin the su- 
perposition of absorptions by gases with different 
pressures and thus only sense an average pressure 
in each case. In ^10. 1.31 we derive a range of possi- 
ble characteristic turbulent Mach numbers for the 
C I-bearing gas 



We explore in £U0.2I some time constants for 
various physical processes that are relevant to our 
work. For instance, the excitation temperatures 
of the two lowest rotational levels of H2 play a 
role many of our pressure determinations. In M 10 . 31 
we find that only on scales of order 100— 1000 AU 
are lag times of any importance in weakening the 
coupling of H2 rotation temperatures to the lo- 
cal kinetic temperature. The same applies to the 
equilibrium between heating and cooling of the 
gas: compressions and decompressions of the gas 
should follow the equilibrium relationship except 
on the smallest scales where the behavior becomes 
more adiabatic in nature. The time scales for 
the equilibria of fine-structure populations and the 
balance between C I and C II are quite short and 
hence apply on all of the relevant size scales. 

In a brief departure from the discussion of tur- 
bulence as a source of pressure deviations, we con- 
sider in EjlO. 41 the possibility that the upper end 
of our pressure distribution function is consistent 
with random interceptions of supernova remnants 
in different stages of development. Finally, we 
summarize our conclusions in fJTT] 

2. Observations 

Our earlier survey (JT01) covered only 21 stars 
that were observed with the guaranteed observ- 
ing time granted to the STIS instrument defi- 
nition team. In order to maximize the observ- 
ing efficiency, the target stars in this study were 
located within two Galactic longitude intervals, 
ones where the HST continuous viewing zoneqj 
intersected the Galactic plane (99° < i < 138°, 
254° < I < 313°). Many additional observations, 
most of which were performed after the study by 
JT01, were more broadly distributed in the sky 
and were taken before the 5-year hiatus of observ- 
ing brought about by STIS instrument failure in 
August 2004. We downloaded from MAST (Multi- 
mission Archive at the Space Telescope Science In- 
stitute) virtually all of the observations performed 
at wavelengths that covered two or more C I tran- 
sitions in the E140H and E230H modes. Once 



3 The continuous viewing zones (CVZs) are two declination 
bands centered on S = ±61? 5 where observations can be 
performed with high efficiency because the targets are not 
occulted by the Earth as the satellite progresses along its 
orbit. 
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again, we made use of the broad wavelength cov- 
erage of STIS to examine as many multiplets as 
possible. These spectra had a resolving power in 
radial velocity equal to 2.6 km s _1 (or 1.5 km s _1 
for the stars observed by JT01 because a nar- 
rower entrance slit was used) (Proffitt et al. 2010). 
Prominent among these observations were those 
performed by a SNAP (snapshot) program con- 
ducted by J. T. Lauroesch (program nr. 8241) in 
1999 and 2000. In the current new study, we have 
also reanalyzed the data presented by JT01 be- 
cause we have now adopted some new, more re- 
fined analysis procedures (see Appendix [A). All of 
the data were processed in the manner described 
by JT01, except that for observations outside their 
survey we did not need to implement an intensity 
rebalancing between MAMA half pixels (see their 
§4.2 for details), since these half pixel intensities 
were binned together beforehand. 

A small fraction of the observations had to be 
rejected because either (1) there was an insuffi- 
cient amount of C I present to perform a meaning- 
ful analysis with the signal-to-noise ratio at hand, 
or (2) the projected rotational velocity of the star 
was so low that stellar features interfered with the 
interstellar ones or made the continua too difficult 
to model. These unsuitable sight lines are iden- 
tified by their target star names in Table [T] We 
also rejected the central stars of planetary nebulae, 
since the interstellar components could be contam- 
inated by contributions from gas in the nebular 
shell. 



Initially, we had considered using observations 
recorded at lower resolution with the E140M mode 
of STIS to broaden the selection of targets, but 
a comparison of the results for a few stars that 
were also observed with the E140H mode indicated 
that unreliable results emerged from the lower res- 
olution data as a result of improper treatments 
of unresolved, saturated profiles. (The correction 
scheme developed by Jenkins (1996) could not be 
used because various lines overlap each other, and 
the optical depths must go through a complicated 
transformation to obtain unique answers for the 
column densities of the three fine-structure levels, 
as described in §5.2.1 of JT01.) 

Table [2] presents the information on the 89 
sightlines included in the present study. They 



Table 1: Rejected Sightlines 
Insufficient Stellar Line 
C I Confusion 



BD+25D2534 CPD-64D481 



BD— 03D2179 


HD 1909 


HD 1999 


HD 3175 


HD 6456 


HD 30122 


HD 6457 


HD 37367 


HD 23873 


HD 43819 


HD 32039 


HD 44743 


HD 64109 


HD 52329 


HD 79931 


HD 62714 


HD 86360 


HD 93237 


HD 92536 


HD 94144 


HD 164340 


HD 106943 


HD 192273 


HD 108610 


HD 195455 


HD 175756 


HD 196867 




HD 201908 




HD 233622 





span path lengths that range from about 0.2 kpc to 
6 kpc and have a median length of 1.9 kpc. We pro- 
cessed all of the data that we felt were acceptable, 
according to the principles outlined in the above 
two paragraphs. We made 2416 separate measure- 
ments, but since we oversampled the wavelength 
resolution of the spectrograph by a factor of 5.3, 
our determinations actually represent only about 
460 independent samples in radial velocity. 

We refrained from applying any special selec- 
tion criteria to make our sampling of regions more 
evenhanded. For this reason, one must be aware of 
certain selection biases in the composite informa- 
tion presented in S}3]below. The following are some 
noteworthy considerations about our sample: 

1. All sightlines terminate at the location of 
a bright, early-type star. Thus, it is in- 
evitable that we will be probing an environ- 
ment near such a star, or in fact a location 
near a grouping of many such stars, since 
they tend to be strongly clustered in space. 
In a number of cases, strong elevations in 
thermal pressures seen in certain radial ve- 
locity channels probably arise from either 
the effects of stellar winds or rapidly expand- 
ing H II regions. As we will show in £19.11 a 
large portion of the C I that we observe re- 



5 



sides in regions with much higher than nor- 
mal densities of starlight radiation. This is 
probably a consequence of the fact that the 
observations are biased toward a sampling 
of the progenitorial cloud complexes out of 
which the stars had formed. 

Regions where the density is low enough (or 
the local temperature or radiation density 
high enough) to shift the ionization equilib- 
rium of carbon atoms more strongly than 
usual to its ionized form are missed in our 
sample. As we show in Figure [I] practically 
all of what is classically known as the warm 
neutral medium (WNM, with T ~ 9000 K) is 
invisible to us; our survey is restricted to the 
phase called the cold neutral medium (CNM, 
with T ~ 80 K), with possibly some very lim- 
ited sensitivity to gas in the thermally unsta- 
ble intermediate temperatures. In addition, 
there are situations where C I is detected at 
certain velocities, but the quality of the data 
is insufficient to measure reliably the ther- 
mal pressures, as we discuss in some detail 
in §A.2I Hence, such regions are excluded. 

Regions of moderate size that are very dense 
will have enough extinction in the UV to 
make stars behind them too faint to ob- 
serve. This effect will result in our miss- 
ing clouds that happen to be strongly com- 
pressed by turbulence or gravity. It is clear 
from the results shown in Column 5 of Ta- 
ble [2] that a cutoff of our sample corresponds 
to a B — V color excess of about 0.5, which 
in turn translates approximately to Af(H) = 
3 x 10 21 cm~ 2 if we use the standard rela- 
tion between E(B-V) and N(R) in the ISM 
(Bohlin et al. 1978; Rachford et al. 2009). 
Also, portions of some of our strongest C I 
absorption profiles were rejected from con- 
sideration because we sensed that they had 
velocity substructures that were saturated 
and not resolved by the instrument. 

Stars in certain programs were selected by 
observers because they had interesting prop- 
erties. Of special relevance to our ther- 
mal pressure outcomes would be the ob- 
servations that were designed to probe re- 
gions that were known to be either disturbed 
(e.g., showing high velocity gas) or at higher 



than normal densities (e.g., showing unusu- 
ally strong molecular absorptions). Often, 
one can sense the characters of such selec- 
tions by reading the abstracts of the pro- 
grams that made the observations^ 



4 The archive root names listed in Column 9 of Table [2] can 
be used as a guide on the MAST HST search web page to 
find the Proposal ID and its abstract. 
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Table 2 
Properties of the Sightlines 



Target 


Galactic Coordinates (deg.) 


Spectral 




uistance 


-T12 J 01 




Archive Exposure 


oil ax 


£ 


b 


lype 


PVR 






ner. 


Root Name(s) 


\ L ) 


(2) 


(3) 


( A\ 




(o) 


{<) 


("J 




CPD-59 2603 


287.590 


-0.687 


05V((f)) 


0.36 


3.5 


77 


6 


O40P01D6Q 
O4QX03010-30 


HD 108 


117.928 


1.250 


06pe 


0.42 


3.8 






O5LH01010-80 


HD 1383 


119.019 


-0.893 


Bill 


0.37 


2.9 






O5C07C010 


HD 3827 


120.788 


-23.226 


B0.7Vn 


0.05 


1.8 






054309010-30 
O54359010-30 


HD 15137 


137.462 


-7.577 


09.5II-IIIn 


0.24 


3.5 


104 


7 


O5LH02010-80 


HD 23478 


160.765 


-17.418 


B3IV 


0.20 


0.47 


55 


7 


O6LJ01020 


HD 24190 


160.389 


-15.184 


B2 Vn 


0.23 


0.82 


66 


7 


O6LJ02020 


HD 24534 (X Per) 


163.083 


-17.137 


09.5III 


0.31 


2.1 


57 


6 


O66P02010 
064813010-20 
O66P01010-20 


HD 27778 (62 Tau) 


172.764 


-17.393 


B3V 


0.34 


0.23 


55 


2 


O59S01010-20 


HD 32040 


196.071 


-22.605 


B9Vn 


0.00 


0.16 






O56L04010-30 
O8MM02010-30 


HD 36408 


188.498 


-8.885 


B7IV 


0.11 


0.19 






O8MM04020-30 


HD 37021 (8 1 Ori B) 


209.007 


-19.384 


B3V 


0.42 


0.56 






O59S02010 


HD 37061 Ori) 


208.926 


-19.274 


B0.5V 


0.44 


0.64 






O59S03010 


HD 37903 


206.853 


-16.538 


B1.5V 


0.29 


0.83 


68 


8 


O59S04010 


HD 40893 


180.086 


4.336 


B0 IV 


0.31 


3.1 


78 


8 


O8NA02010-20 


HD 43818 (11 Gem) 


188.489 


3.874 


B0 II 


0.45 


1.9 






O5C07I010 


HD 44173 


199.002 


-1.316 


B5III 


0.05 


0.52 






O5C020010 


HD 52266 


219.133 


-0.680 


09.5 IVn 


0.22 


1.8 






O5C027010 


HD 69106 


254.519 


-1.331 


B0.5IVnn 


0.14 


1.5 


80 


10 


O5LH03010-50 


HD 71634 


273.326 


-11.524 


B7IV 


0.09 


0.32 






O5C090010 


HD 72754 (FY Vel) 


266.828 


-5.815 


B2I:pe 


0.31 


3.9 






O5C03E010 


HD 75309 


265.857 


-1.900 


BlIIp 


0.18 


2.9 


65 


3 


O5C05B010 


HD 79186 (GX Vel) 


267.366 


2.252 


B5Ia 


0.23 


1.9 






O5C092010 


HD 88115 


285.317 


-5.530 


B1.5IIn 


0.12 


3.7 


145 


3 


054305010-60 


HD 91824 


285.698 


0.067 


07 V 


0.22 


3.0 


61 


6 


O5C095010 


HD 91983 


285.877 


0.053 


B1III 


0.14 


3.0 


61 


7 


O5C08N010 


HD 93205 


287.568 


-0.706 


03Vf+ 


0.34 


3.3 


105 


6 


O4QX01010-40 



Table 2 — Continued 



Target 


Galactic Coordinates (deg.) 


Spectral 




Distance a 


H2 Toi b 




Archive Exposure 


Star 


o 
I 





Type 


E(B - V) a 


(kpc) 


(K) 


Ref. c 


Root Name(s) 


(1) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HD 93222 


287.738 


-1.016 


07IIK 


0.32 


3.6 


77 


6 


O4QX02010-40 


HD 93843 


288.243 


-0.902 


05 Hlf 


0.24 


3.5 


107 


10 


O5LH04010— 40 


HD 94454 


295.693 


-14.725 


B8III 


0.19 


0.30 


74 


7 


O6LJ0H010 


HD 94493 


289.016 


-1.177 


Bl lb 


0.15 


3.4 






054306010—20 


HD 99857 


294.779 


-4.940 


B0. 5 lb 


0.27 


3.5 


83 


10 


054301010-60 


















054301020 


















054301030 


















054301040 


















054301050 


















054301060 


HD 99872 


296.692 


-10.617 


B3 V 


0.29 


0.24 


66 


7 


O6LJ0I020 


HD 102065 


300.027 


-17.996 


B2 V 


0.28 


0.18 


59 


6 


O4O001010-30 


HD 103779 


296.848 


-1.023 


B0. 5 lab 


0.17 


4.3 


86 


6 


054302010—20 


HD 104705 (DF Cru) 


297.456 


-0.336 


B0 lb 


0.17 


5.0 


92 


6 


O57R01010, 30 


HD 106343 (DL Cru) 


298.933 


-1.825 


B1.5Ia 


0.23 


3.3 






054310010—20 


HD 108002 


300.158 


-2.482 


B2Ia/Iab 


0.18 


4.2 


77 


7 


O6LJ08020 


HD 108639 


300.218 


1.950 


B0. 2 III 


0.26 


2.4 


88 


7 


O6LJ0A020 


HD 109399 


301.716 


-9.883 


B0. 7 II 


0.19 


2.9 






054303010-20 


HD 111934 (BU Cru) 


303.204 


2.514 


B1.5 lb 


0.32 


2.3 






O5C03N010 


HD 112999 


304.176 


2.176 


B6 Vn 


0.17 


0.45 


96 


7 


O6LJ0C010— 20 


HD 114886 


305.522 


-0.826 


09IIIn 


0.32 


1.8 


92 


7 


O6LJ0D020 


HD 115071 


305.766 


0.153 


B0. 5 Vn 


0.40 


2.7 


71 


7 


O6LJ0E010— 20 


HD 115455 


306.063 


0.216 


07.5 III 


0.40 


2.6 


81 


7 


O6LJ0F010-20 


HD 116781 


307.053 


-0.065 


B0 Bine 


0.31 


2.2 






O5LH05010-40 


HD 116852 


304.884 


-16.131 


09 III 


0.14 


4.5 


70 


6 


O5C01C010 


















063571010 


















O8NA03010-20 


HD 120086 


329.611 


57.505 


B2V 


0.04 


0.99 






O5LH06010-50 


HD 121968 


333.976 


55.840 


Bl V 


0.11 


3.1 


38 


6 


O57R02010-20 


HD 122879 


312.264 


1.791 


B0 la 


0.29 


3.3 


90 


7 


O5C037010 



O5LH07010-40 
O6LZ57010 



Table 2 
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Target 


Galactic Coordinates (deg.) 


Spectral 




Distance a 


H2 Toi b 




Archive Exposure 


Star 


n 
t 


r 




Type 


E(B - Vf 


(kpc) 


(K) 


Ref. c 


Root Name(s) 


(1) 


( 2 ) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HD 124314 


312.667 


-0.425 


06 Vnf 


0.43 


1.4 


74 


7 


054307010-20 


HD 140037 


340.151 


18.042 


B5 III 


0.08 


0.77 






O6LJ04010 


HD 142315 


348.981 


23.300 


B9V 


0.10 


0.15 






O5C03Y010 


HD 142763 


31.616 


46.960 


B8III 


0.01 


0.28 






O5C040010 


HD 144965 


339.044 


8.418 


B2 Vne 


0.27 


0.51 


70 


7 


O6LJ05010 


HD 147683 


344.858 


10.089 


B4 V+B4 V 


0.28 


0.37 


58 


7 


O6LJ06020 


HD 147888 (p Opn D) 


353.648 


17.710 


B3 V 


0.42 


0.12 


44 


7 


O59S05010 


HD 148594 


350.930 


13.940 


B8 Vnn 


0.18 


0.19 






O5C04A010 


HD 148937 


336.368 


-0.218 


06.51 


0.55 


2.2 






O6F301010-20 


HD 152590 


344.842 


1.830 


07 V 


0.37 


3.6 


64 


7 


O5C08P010 
O8NA04010— 20 


HD 156110 


70.996 


35.713 


B3 Vn 


0.03 


0.62 






O5C01K010 


HD 157857 


12.972 


13.311 


06.5 IHf 


0.37 


3.1 


86 


7 


O5C04D010 


HD 165246 


6.400 


-1.562 


08 Vn 


0.33 


1.9 






O8NA05010— 20 


HD 175360 


12.531 


-11.289 


B6III 


0.12 


0.24 






O5C047010 


HD 177989 


17.814 


-11.881 


B0 III 


0.11 


6.0 


52 


6 


O57R04010-20 
O57R03010-20 


HD 185418 


53.604 


-2.171 


B0. 5 V 


0.38 


1.2 


105 


6 


O5C01Q010 


TTT\ -1 AA/>OA 

HD 192639 


74.903 


1.479 


07Ibf 


0.56 


2.1 


98 


2 


O5C08T010 


HD 195965 


85.707 


4.995 


B0 V 


0.19 


1.1 


91 


7 


O6BG01010— 20 


HD 198478 (55 Cyg) 


85.755 


1.490 


B3 la 


0.43 


1.3 






O5C06J010 


HD 198781 


99.946 


12.614 


B0. 5 V 


0.26 


0.69 


65 


7 


O5C049010 


HD 201345 


78.438 


—9.544 


09 V 


0.14 


2.2 


147 


6 


O5C050010 
O6359P010 


HD 202347 


88.225 


-2.077 


B1.5 V 


0.11 


0.95 


116 


3 


O5G301010, 40-50 


HD 203374 


100.514 


8.622 


B2Vn 


0.43 


0.34 


87 


10 


O5LH08010-60 


HD 203532 


309.461 


-31.739 


B3IV 


0.24 


0.22 


49 


6 


O5C01S010 


HD 206267 


99.292 


3.738 


06.5V 


0.45 


0.86 


65 


2 


O5LH09010-40 


HD 206773 


99.802 


3.620 


B0V:nnep 


0.39 


0.82 


94 


4 


O5C04T010 


HD 207198 


103.138 


6.995 


09.5Ib-II 


0.47 


1.3 


66 


2 


O59S06010-20 


HD 208440 


104.031 


6.439 


B1V 


0.27 


1.1 


75 


4 


O5C06M010 


HD 208947 


106.550 


8.996 


B2V 


0.16 


0.56 






O5LH0A010-40 



Table 2 — Continued 



Target 


Galactic Coordinates (deg.) 


Y-S v_- \j L L CX>± 




Distance a 


H2 Tq\ 




Archive Exposure 


Star 




u 




Tvnp 


e(b — vy 




(K) 


Ref. c 


Root TVflTTIpf ^ 1 




V 2 ) 


(3J 


(4) 


\°) 




(7) 
\ < ) 


(H) 




HD 209339 


104.579 


5.869 


BO IV 


0.24 


1.2 


90 


4 


O5LH0B010-40 


















O6LZ92010 


HD 210809 


99.849 


-3.130 


09Iab 


0.28 


4.3 


87 


7 


O5C01V010 


HD 210839 (A Cep) 


103.829 


2.611 


06Infp 


0.49 


1.1 


72 


2 


054304010-20 


HD 212791 


101.644 


-4.303 


B3ne 


0.18 


0.62 






O5C04Q010 


HD 218915 


108.064 


-6.893 


09.5 Iabe 


0.21 


5.0 


86 


6 


O57R05010, 30 


HD 219188 


83.031 


-50.172 


B0.5IIIn 


0.09 


2.1 


103 


1 


O6E701010 


















O8DP01010 


















O8SW01010 


HD 220057 


112.131 


0.210 


B3IV 


0.17 


0.77 


65 


7 


O5C01X010 


HD 224151 


115.438 


-4.644 


B0.5II-III 


0.34 


1.3 


252 


10 


054308010-20 


HDE 232522 


130.701 


-6.715 


Bill 


0.14 


6.1 






O5C08J010 


HDE 303308 


287.595 


-0.613 


03 Vf 


0.33 


3.8 


86 


6 


O4QX04010-40 



a B — V color excesses and distances to the stars were either taken from listings of the same stars in Bowen et al (2008) or Jenkins 
(2009), or else they were computed by using the same procedures that they invoked. 

The molecular hydrogen rotational temperature from J = to 1 that was adopted as an indicator for the kinetic temperature 
of the intervening gas. 

c Reference for the source of the T01 value given in the previous column: (1) Savage et al. (1977); (2) Rachford et al. (2002); (3) 
Andre et al. (2003); (4) Pan et al. (2005); (5) Lee et al. (2007); (6) Burgh et al. (2007); (7) Sheffer et al. (2008); (8) Rachford et 
al. (2009); (9) Burgh et al. (2010); (10) "J. M. Shull (2009) in preparation" listed in Burgh et al. (2010); (11) Jensen et al. (2010). 




-2-10123 

Log n(H) 

Fig. 1. — A plot of the thermal pressure, log(p/fc), vs. the density of hydrogen nuclei, log n(H), showing 
the locations where heating equals cooling in the ISM near our part of the Galaxy (thick curve), according 
to the thermal equilibrium calculations by Wolfirc ct al. (2003) (their "standard model"; see their Fig. 8). 
Portions of this curve where the slopes are positive are thermally stable and form the distinct phases called 
the warm neutral medium (WNM) and cold neutral medium (CNM), as indicated. The portion of the curve 
that has a negative slope has a balance between heating and cooling, but is thermally unstable (Field 1965). 
In the absence of rapidly changing pressures and densities due to turbulence, the lowest allowable pressure 
for the CNM is at the horizontal dash dot line labeled "CNM log(p//c) m i n ." Different temperatures in this 
diagram are revealed by the straight, dashed lines, constructed using the assumption that He/H=0.09 and, 
for T < 1000 K, /(H2) = 0.6 (see H4.3[) . The thin, gently curved lines show constant values for the expected 
values of C Itotai/(C II + C Itotai), as indicated, according to our equation for ionization equilibrium (see 
Eq. [T]in ^4.41 and the accompanying text), under the assumption that the starlight intensity is equal to the 
average level Jo given by Mathis et al. (1983). These curves demonstrate that the WNM is virtually invisible 
in our survey of C I. 
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3. CI Results 

Descriptions of our analysis and the mathemat- 
ical details of the interpretation of the C I absorp- 
tion multiplcts covered by HST were presented by 
JT01. Except for some enhancements in technique 
and the use of more up to date atomic data dis- 
cussed in Appendix El we have implemented once 
again the methods of JT01. 

Briefly, after we normalized the intensity pro- 
files to an assumed continuum (that usually varies 
smoothly with wavelength) , we converted them to 
apparent optical depths (Savage & Sembach 1991). 
For the average spread in radial velocity of the 
C I in a typical sightline, the individual lines in 
any given multiplet overlap each other. This intro- 
duces confusion in the interpretation of the optical 
depths. However, we can unravel this confusion 
by observing different multiplets, because the lo- 
cations of different transitions with respect to each 
other change, thus allowing one to resolve ambigu- 
ous mixtures of opacities. JT01 devised a way to 
construct a system of linear equations that could 
be solved to reveal the apparent column densities 
N(C I), N(C P), and N(C I**) as a function of 
velocity!! Once this has been done, we can eval- 
uate the quantities /l = N(C l*)/N(C Itotai) and 
/2 = N(C l**)/N(C hotai) , which are useful rep- 
resentations of the excitation conditions when we 
want to understand not only the physical condi- 
tions in any given absorbing region, but also pos- 
sible combinations of contributions from differing 
regions that overlap each other at a particular ve- 
locity. 

As explained originally by Jenkins & Shaya 
(1979) and once again by JT01, the balance of 
collisional excitations (and de-excitations) against 



5 The notation adopted here is consistent with that of JT01: 
N(C I) refers to the column density of atomic carbon in 
its 3 Po ground fine-structure state, while N(C I*) and 
JV(C I**) refer to the column densities of the excited 3 Pi 
and 3 P2 levels, respectively. The quantity 7V(C Itotai) 
equals the sum of the column densities in all three levels. 
Strictly speaking, we measure apparent column densities 
[N a in the notation of Savage & Sembach (1991)], which 
differ from true column densities because the recorded in- 
tensities are smoothed by the instrumental line spread func- 
tion. In the interest of simplicity, we will refer to such ap- 
parent column densities as simply N and treat them is if 
they were true column densities. Possible errors in this as- 
sumption and our avoidance of cases where they are large 
are discussed in EIA.2.41 



the spontaneous radiative decay of the excited lev- 
els establishes an equilibrium value for the level 
populations that depends on the local density and 
temperature (and to a much lesser extent, the 
composition of the gas). As the densities increase 
at any given temperature, the locations of points 
on a diagram of fl vs. /2 trace an upward arching 
curve (see Fig. [2]) that stretches from the origin fol- 
low densities to a point at very high densities that 
approaches a Boltzmann distribution for the levels 
at the temperature in question. 

When the absorptions from two or more regions 
are superposed at a single velocity, the outcome 
for fl and /2 is at the "center of mass" for the 
values that apply to the individual contributors, 
with respective weights equal to their values of 
N(C Itotai)- This outcome is not without some 
ambiguity, since various combinations of condi- 
tions in any ensemble of different clouds can pro- 
duce the same result. We will address this issue 
later in Sj5]when we make a simplifying assumption 
about such mixtures, and in i)9.3l we will discuss 
the consequences of possible averaging effects that 
are difficult to recognize. Some additional com- 
plexity emerges when one considers the effects of 
optical pumping by starlight photons, which we 
will cover in $6] 

Figure [5] shows the outcome for all of our mea- 
surements of (fl, /2) at each velocity interval that 
showed acceptable results, and for every star in the 
survey. The area of each dot in this diagram is pro- 
portional to N(C Total)- It is clear that practically 
all of the points fall above the equilibrium calcula- 
tions for fl and /2, even for a 300 K temperature 
that is well above the nominal values for the CNM. 
From this we conclude that either there is always 
a mixture of two or more regions with vastly dif- 
fering conditions for every velocity channel or that 
for some reason(s) the curves are incorrect or in- 
appropriate (we will touch upon this issue later in 
£110.5.11) . A generalized picture of how we inter- 
pret some plausible admixtures will be presented 
in 
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Fig. 2. — Measurements of fl and /2 for all velocity bins, each having a width of 0.5 km s~ , for the 89 stars 
in the survey that had uncertainties and cr(/2) less than 0.03. The area of each dot is proportional 

to the respective value of N(C Itotai), with a normalization in size as shown in the box at the upper left 
portion of the plot. The white x located at fl = 0.209, /2 = 0.068 represents the "center of mass" of all 
of the dots. The curves indicate the expected level populations for three different temperatures, assuming 
the gas mixture as specified in ^4.31 with different values of \og(p/k) indicated with dots (adjacent dots 
represent differences of 0.1 dex). The large open circles on the curves indicate integer values of \og(p/k) 
with accompanying numbers to indicate their values. Populations that are proportional to the degeneracies 
of the levels are indicated by the + sign labeled "(p, T) — »• oo." 
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4. Properties of the C I-Bearing Gas 

4.1. Corrections for the Ionization of Car- 
bon 

It is important to realize that C I in the gas 
that we observe is a minor constituent, since the 
ionization potential of neutral carbon (11.26 eV) is 
below that of neutral hydrogen (13.6 eV). Most of 
the carbon atoms in H I regions are singly ionized. 
The relative proportion of carbon in the neutral 
form can vary by enormous factors according to 
the local density of the gas and the strength of 
the local radiation field that is responsible for ion- 
izing the atoms. As we attempt to derive an even- 
handed picture of the pressure distribution for the 
general neutral gas, rather than just a value that is 
weighted in proportion to the amount of C I that 
is present, we must devise a means for asssessing 
how much C II accompanies the C I. In essence, 
we use C II as an indicator for the total amount 
of the neutral material. 

Direct measures of N(C II) are very difficult 
to carry out. The only available transition in the 
wavelength bands covered by our survey is the one 
at 1334.53 A. This line is strongly saturated, and 
the only way to measure N(C II) with this feature 
is by sensing the strength of its damping wings 
(Sofia et al. 2011). However, this measure applies 
to C II at all velocities, rather than at the veloci- 
ties where we are able to make use of information 
from C I. While there exists a very weak intersys- 
tem transition at 2325.40 A, this feature is outside 
the wavelength coverage of most of our observa- 
tions and also requires a very high signal-to-noise 
ratio for a reliable detection (Sofia et al. 2004). 

To overcome our inability to measure directly 
N(C II) as a function of velocity, we instead used 
O I as a proxy for C II. The very weak O I in- 
tersystem line at 1356 A is ideal for tracing all 
but the smallest column densities of material per 
unit velocity. For velocity intervals over which the 
C I absorptions could be measured reliably, we 
found that only on very rare occasions was the 
O I line too weak to measure. For such instances 
we had to use the weakest line of S II at 1250 A 
as a substitute for O 10 For deriving N(C II), 



' Over very restricted velocity intervals there was an inter- 
mediate range of column densities per unit velocity where 
the O I line was too weak to observe (apparent optical 



we assumed that C, O and S were depleted below 
their respective protosolar abundances (Lodders 
2003) by amounts equal to —0.162 dex, —0.123 dex 
and —0.275 dex, respectively, which corresponds 
to a moderate depiction strength (F* = 0.5) in 
the generalized representation of Jenkins (2009). 
The assumed abundance of O relative to C could 
be in error by about 0.05 dex if actual the deple- 
tion strength is either = 0.0 or 1.0 instead of 
0.5. Also, a few new determinations of iV(C II) 
by Sofia et al. (2011) based on fitting the damp- 
ing wings of the strong line at 1334.53 A instead of 
using the weak intersystem line suggest that the 
abundances of C in the ISM may be about 0.3 dex 
lower than stated above, but this change would 
create a uniform offset that would apply to all of 
our cases. 

One can imagine that some of the C Il-bearing 
gas may be situated in fully ionized regions inter- 
sected by our sight lines, especially if most of the 
ionizing radiation in the ISM is below the ioniza- 
tion potential of singly ionized carbon (24.38 eV). 
At velocities where we rely on O I as a proxy, we 
are confident that our estimate for the amount of 
C II applies only to neutral gas, since H II regions 
are devoid of O I because the ionizations of O and 
H are strongly coupled by a charge exchange reac- 
tion with a large rate constant (Field & Steigman 
1971; Chambaud et al. 1980; Standi et al. 1999). 
The same is not true for S II; its behavior should 
be similar to that of C II (singly ionized S and 
C have ionization potentials within 1 eV of each 
other). However, we find that the continuity over 
velocity between the O I and S II profiles is usually 
good, which argues against the existence of much 
contamination from H II regions. 

4.2. Kinetic Temperature of the C I- 
bearing Gas 

The relationship between the thermal pressure 
and the fine-structure excitations has a weak de- 
pendence on temperature. Thus, even though this 
effect is small, we must still try to reduce the ambi- 
guities in the measurement of p when either n(H) 
or T is unknown. Fortunately, we can make use 



depth T a < 0.05) and the S II line was badly saturated 
(r a > 2.5). In this range, we adopted a geometric mean as 
an approximate compromise between the respective upper 
and lower limits. 
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of the fact that molecular hydrogen usually ac- 
companies C I and thus we can utilize measure- 
ments of the J = to 1 rotation temperature of 
H2 to indicate the most probable kinetic tempera- 
ture of the material for which we are making pres- 
sure measurements. Columns (7) and (8) of Ta- 
ble [2] show the outcomes of such rotation temper- 
ature measurements, Tqi, along with the sources 
in the literature where such measurements were 
listed. In cases where Tqi was not measured, we 
had no alternative but to adopt an arbitrary num- 
ber T = 80 K, which is close to the median value 
of all of the measurements. However, variations in 
Tqi from one sightline to the next have rms devi- 
ations of only 30 K, and deviations of this magni- 
tude alter the outcome for log(p/k) by only about 
0.06 dex. 

We acknowledge the presence of an unavoidable 
limitation that Tqi shows simply an average over 
all velocities, and it is weighted in proportion to 
the local density of hydrogen molecules, n(H2), in- 
stead of n(C Itotai)- In some circumstances varia- 
tions in temperature across different regions along 
our sightlines might compromise the accuracy of 
our results, but this is probably not a very impor- 
tant effect. 

4.3. The Particle Mix 

The composition of the gas has a small, but 
nonnegligible effect on the expected outcomes for 
/l and /2. For instance, for 2.8 < log(p/fc) < 3.8 
the inferred pressure for a given (/l, /2) for pure 
atomic hydrogen is about 0.1 dex lower than for 
the equivalent number density of pure H2 with 
Tqi = 80 K. Thus, in order to minimize the error in 
the interpretation of the measurements, it is good 
to adopt an estimate for the most probable mix of 
gas constituents. Any deviations in the true con- 
ditions from whatever we adopt for the molecular 
fraction, /(H 2 ) = 2n(H 2 )/[2n(H 2 ) + n(H I)], will 
result in an error for \og(p/k) of less than 0.1 dex. 

We estimate that approximately half of the ma- 
terial in our lines of sight arise from the WNM, 
which is free of H2, while the remaining half 
(CNM) that we can see with C I has an appre- 
ciable molecular content. If we assume that the 
CNM has /(H 2 ) = 0.60, then an overall value 
/(H 2 ) = 0.42 would apply to the entire sightline. 
The latter value is very close to the median out- 
come for J(H 2 ) found by Rachford et al. (2009) 



in their FUSE survey of sightlines similar to the 
ones in the present study. We therefore adopt the 
assumption that /(H2) = 0.60 in our C I-bearing 
gas, and the mix of ortho- and para-H2 is gov- 
erned by the determination of X01 (which is set to 
80 K if unknown). Another constituent is helium, 
whose ratio to hydrogen in atoms and molecules is 
assumed to be the protosolar value 0.094 given by 
Lodders (2003). In their normal concentrations in 
the CNM, electrons and protons have a negligible 
influence on fl and /2. 

4.4. Local Starlight Intensity 

In H4.ll we explained how we derive the amount 
of C II that accompanies the C I. This determi- 
nation at any particular velocity in a given sight- 
line has two applications. First, as we indicated 
earlier, it represents our best estimate for the to- 
tal amount of neutral material associated with the 
C I. Second, the ratio of C II to C I allows us to 
estimate the local radiation density, which in turn 
will be useful for applying a correction to the equi- 
librium fl and /2 values that allows for effects of 
optical pumping of the fine-structure levels (see 
H6.ll and ijA.41 for details). As we will show later 
in i)9.1l and ^10.5.41 the radiation density outcomes 
are by themselves of special interest in our overall 
outlook on the distribution of thermal pressures. 

The ionization balance of carbon atoms in the 
neutral ISM is given by the relation 

n(C I tot al)(I//o)r (C I) = n(C II)[a e (C ll,T)n(e) 

+a g (C n,n(e),/,T)n(H)] (1) 

where the photoionization rate To(C I) = 2.0 x 
10~ 10 s _1 (Weingartner & Draine 2001a) if the ra- 
diation field density / is equal to a value Iq spec- 
ified by Mathis et al. (1983) for the average in- 
tensity of ultraviolet starlight in our part of the 
Galaxy, a e (C II, T) is the radiative plus dielec- 
tronic recombination coefficient of C II with free 
electrons as a function of temperature T (Shull 
& Van Steenberg 1982), and a g (G II, n(e), I, T) 
is the C II recombination rate due to collisions 
with dust grains (and subsequent transfer of an 
electron) normalized to the local hydrogen den- 
sity (Weingartner & Draine 2001a) To solve for 



7 In the notation of Weingartner & Draine (2001a), this elec- 
tron transfer rate from grains is expressed as a g (C+, i/>, T), 
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n(e), we assume that free electrons are created by 
both the photoionization of some heavy elements, 
amounting to 2 x 10~ 4 n(ff ), supplemented by elec- 
trons liberated from the cosmic-ray ionization of 
hydrogen at a rate Ccr = 2x 10~ 16 s -1 (Indriolo 
ct al. 2007; Neufeld et al. 2010) @ In order to 
estimate the density of electrons created by these 
cosmic-ray ionizations, one must calculate the bal- 
ance between the creation of free protons with a 
density n(p) against their recombination with free 
electrons and also electrons on dust grains using 
an equation analogous to Eq. [TJ 

CcRn(H)=n(p)[a e (H U,T)n(e) 
+a g (RU,n(e),I,T)n(K)} . (2) 

As with the case for C II, we obtain a formula 
for a s (H II, n(e), I, T) from Weingartner & Draine 
(2001a). 

There may be some shortcomings in our simple 
formulation in Eq. Q] Welty et al. (1999) found 
inconsistencies in the determinations of electron 
densities using the ratios of neutral and ionized 
forms of different elements, and these problems are 
not resolved when the grain recombination pro- 
cesses are included (Weingartner & Draine 2001a). 
Either the rates incorporated into Eq. Q] are inac- 
curate or other kinds of reactions may be impor- 
tant, such as charge exchange with protons, or the 
formation and destruction of CO (van Dishoeck & 
Black 1988) or other C-bearing diatomic molecules 
(Prasad & Huntress 1980; van Dishoeck & Black 
1986, 1989). 

If we use our determination of N(C II) dis- 
cussed in ^4.11 compare it with iV(C Itotai), and 
apply the equilibrium condition expressed by Eq.Q] 
to determine the local starlight density, we should 
obtain a reasonably accurate result for I/Iq pro- 
vided that there is not a large amount of C II 
arising from the WNM at exactly the same veloc- 
ity. If such an additional contribution is present, 
we will overestimate the radiation density. Like- 
wise, this overestimate of N(C II) accompanying 
the C I will give a disproportionately large figure 
for the total amount of neutral material for the 



where if) = GT2n(e) 1 and G = 1.13 for the interstellar 
radiation field of Mathis et al. (1983). 

For the column densities of hydrogen considered here, the 
average ionization from x-rays is almost negligible by com- 
parison (Wolfire et al. 1995). 



particular measurement at hand. 

5. Admixtures of Different Kinds of Gas 

It is clear that the 2-dimensional distribution 
of the points shown in Fig. [^represents a complex 
mixture of high and low pressure gas. Some in- 
sight on the possible nature of this mixture may 
be gained by stripping away the information about 
the 2-dimensional scatter of the points and focus- 
ing on just a composite value of (fl, /2) for all 
of the measurements, that is, a single "center of 
mass" of all of the points in the diagram. The lo- 
cation of this value is shown by a white x in the 
figure. The question we now ask is whether or not 
some simple, generalized distribution function for 
pressures can reproduce the observed composite 
{fl, /2) pair. 

A very elementary but plausible pressure rela- 
tionship to propose is a lognormal distribution, 
which is appropriate for situations where ran- 
dom pressure fluctuations arise from turbulence 
(Vazquez Semadeni 1994; Nordlund & Padoan 
1999; Kritsuk et al. 2007) - see g9J] Panels 
(a) through (c) in Figure [3] show reconstructions 
of the combinations of fl and /2 for such a dis- 
tribution for three different values for the width 
of the distribution in log(p/fc) and a single value 
for the location of the peak. One might initially 
suppose that the curvature of the high pressure 
tail in this distribution along an arc that traces 
the single-region (fl, /2) combinations could pull 
the distribution's composite (fl, /2) to a loca- 
tion above the curve at low pressures. However, 
Figure [3] illustrates that this curvature appears to 
be insufficient to make our model lognormal com- 
posite values rise as high (in /2) as the measured 
ones without passing beyond the composite mea- 
surement of fl. Instead, it appears that a more 
complex picture is called for, one that requires the 
use of a bimodial distribution of pressures. The 
simplest such model is to propose the existence of a 
separate, small contribution from material at pres- 
sures log(p/fc) > 5.5 and a temperature T > 80 K. 
Panel (d) of the figure shows that such a contribu- 
tion can solve our problem with the elevated com- 
posite (fl, /2). This high pressure contribution 
seems to be present in a majority of the sightlines, 
since most of the individual points that are shown 
in Fig. [2] are pulled above the curves. 
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Fig. 3. — A schematic demonstration of the behavior of the composite values of f\ and /2 when the pressure 
distribution follows a simple lognormal behavior. The distribution is approximated by a discrete collection 
of H I packets spaced 0.1 dex apart in pressure, illustrated by black dots strung along the equilibrium curve, 
with the area of each dot indicating the amount of H I. After factoring in the ionization equilibrium equation 
for carbon atoms, the amounts of C I are strongly biased in favor of higher densities. The amounts of C I 
are indicated by the areas enclosed by open circles. The "center of mass" for the C I packets appears at 
the location of the square with an arrow pointing toward it, while the observed composite (/l, /2) shown 
in Fig. [5] is indicated with an x. Panels (a) through (c) show lognormal distributions with 3 successively 
increasing values of cr, while panel (rf) shows the outcome when a small, additional contribution of very high 
density gas is present (very small circles at the top of the equilibrium curve). 
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Fig. 4. — A demonstration of how an observed (/l, /2) at a given velocity is decomposed into a superposition 
of low and high pressure regions. As in Fig. [2l the equilibrium track is marked with a scale in log(p/fc), 
with open circles and accompanying numbers showing integer values of this quantity. An assigned location 
of (0.38, 0.49) in this diagram applies to the high pressure component, and a line that projects from this 
point through the observed (/l, /2) intersects the equilibrium curve (for a given temperature) at a point x 
that corresponds to the pressure of the low pressure component (in this case \og(p/k)\ ow = 3.5). The relative 
distances along the projection line indicate fractions of C I in the two components: in this depiction, the 
length of the segment above "Obs." indicates that fraction of gas at normal (low) pressures is gi ow = 0.90, 
and the remaining fraction in the high pressure component is indicated by the length of the lower line 
segment, yielding ghigh = 0.10. 
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The exact properties of the distribution of high 
pressure material is not known, but for (/l, /2) 
outcomes not far above the lower portions of the 
equilibrium curves, such details do not matter 
much for the low pressure gas. We need only to 
know the general vicinity in the upper portion of 
the diagram where such material resides. Hence- 
forth, we adopt (0.38, 0.49) for a fiducial fl and 
/2 for the high pressure component (which cor- 
responds approximately to T = 300 K, n(R I) = 
4000 cm -3 ), and consider that from one observa- 
tion to the next, the relative proportion of this 
gas is some small fraction of the total that can 
vary from one case to the next. The consequences 
of the high-pressure reference point differing from 
reality will be discussed briefly in $7] 

Figure U is a schematic illustration of how we 
geometrically decompose an observed combination 
of fl and /2 into a superposition of the two pro- 
posed types of gas, one at a very high pressure 
(but with poorly known physical conditions) and 
the other at a normal, low pressure. The basic 
strategy is to find where a projection from the 
assumed high density locus (0.38, 0.49) through 
the observed combination of fl and /2 extends 
to a specific point on the equilibrium curve drawn 
for an appropriate temperature, as defined by Tqi, 
if available. We regard the pressure that corre- 
sponds to this intersection point to represent the 
proper result for the low density gas. The ratio 
of C I in the high pressure gas to the overall total 
is given by the quantity 1 — g\ ow (or simply ghigh) 
shown in the diagram. 

Two different considerations governed our 
choice of the fiducial high pressure (fl, f2) to 
be situated near the top of the points shown in 
Fig. [2j One is based on the plausibility that for 
most of the ordinary lines of sight the amount of 
this gas (in terms of the total gas, not just C I) 
is probably a small fraction of all of the gas. At 
the highest pressures, C I is more conspicuous 
(because of the shift in the ionization balance to- 
ward the neutral form of C), which in turn leads 
to a smaller quantity for the inferred amount of 
singly-ionized carbon. The other consideration is 
a practical one: the numerical results for the de- 
compositions into low and high pressure gas are 
more stable when the high pressure point is well 
removed from the measured values of fl and /2 
- this becomes important for outcomes at moder- 



ately high pressures above the main distribution 
of low pressures. 

In the next several sections, we will concentrate 
on the measurements of the low pressure com- 
ponent, using the method just described. Later, 
in £110. 5[ we will turn our attention to the small 
amount of gas at high pressures and discuss its 
possible significance in our understanding of pro- 
cesses in the diffuse ISM. 

6. Derivations of Thermal Pressures 
6.1. Initial Estimates of Conditions 

For each measurement of (fl, /2), we apply the 
construction demonstrated in Fig. 3] to determine 
the quantities log(p/fc)i ow and gi ow . This determi- 
nation is based on the initial assumption that the 
radiation field intensity I(X) in the gas in the low 
pressure regime is equal to the average Galactic 
value Io(X) (see £|4.4p . Any deviation of the true 
intensity from this relationship will result in an er- 
ror in \og(p/k) because an incorrect calculation of 
the optical pumping rate was applied. (Different 
pumping rates have virtually no effect on the value 
of ffiow since the equilibrium values of fl and /2 
simply shift along the curve that represents differ- 
ent pressures.) 

In order to obtain a more accurate result, we 
must evaluate how much the true intensity I(X) 
differs from Iq(X), so that our pumping correc- 
tion will be more accurate. We make the sim- 
plifying assumption that the relative distribution 
of intensity over A does not change appreciably 
from one location to the next, but that the overall 
level of radiation is a free parameter that can vary. 
We then estimate an initial approximation for the 
value of this parameter, which we call I/Iq. To 
derive this estimate we make use of the ionization 
equilibrium equation, Eq. [TJ to solve for I/I for 
gas at any particular velocity increment, after re- 
placing n(C II) by N(C II) (as determined accord- 
ing to the method described in i)4.ip . n(C Itotai) 
by gi ow N(C Itotai) and n(H) by our initial approx- 
imate value of p/(kToi). In making the substitu- 
tion for 7i (C Itotai), we assume that it is safe to 
declare that virtually all of the C II is associated 
with the low pressure component of C I, but this 
condition could be violated if the radiation den- 
sity experienced by the high pressure gas is many 
orders of magnitude higher than that of the gas at 
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ordinary pressures. 

6.2. Convergence to Final Values 

After evaluating the new radiation intensity 
level I /la, we are in a position to repeat the cal- 
culation of log(p/fc)iow using a better representa- 
tion for the shifts in the expected values of fl 
and fl caused by optical pumping!^ However, 
the new value for the pressure will have an im- 
pact on the density n(H) used in the equation for 
the ionization equilibrium, hence the calculation of 
this balance must be repeated in order to derive a 
modified number for radiation enhancement factor 
I / Iq , one that is better suited for a more accurate 
derivation of the pressure. We cycle through the 
alternation between pressure and ionization calcu- 
lations many times until the densities and intensi- 
ties converge to stable solutions. 

7. Overview of Sightlines 

Table [3] presents a number of properties of the 
C I and (inferred) C II data that we obtained for 
the sightlines that were suitable for study. The 
numbers in this table give general indications in- 
tegrated over velocity; they were not used in the 
analysis of the pressure distribution, which relied 
on the more detailed results that we obtained for 
the explicit velocity channels. 

We show our estimates of the total column den- 
sities of C II in Column (3) of the table. They com- 
pare favorably with the few direct determinations 
reported in the literature (see note /); rms devi- 
ations between our values and others amount to 
0.22 dex. When we consider that the direct mea- 
surements of N(C II) have quoted errors of order 
0.1 dex, the magnitudes of the disagreements in- 
dicate that our values are probably uncertain by 
about 0.20 dex. The largest deviations seem to oc- 
cur for cases where the other determinations are 
higher than ours, which may indicate that we are 



9 The shift in the outcome for log(p/fc) depends not only 
on the strength of the pumping field intensity, but also 
on p/k itself. Figure 6 of Jenkins & Shaya (1979) shows 
the (/l, /2) equilibrium tracks for / /Iq = 1 and 10, but 
in terms of our revision of the pumping rates derived in 
3A.4I these tracks are equivalent to present-day values for 
I/Iq equal to about 1.5 and 15. Representative values 
for log(p/fc) and log(//7o) are listed for each sight line in 
Columns (6) and (7) of Table [3] 



not registering some C II in ionized gas because 
we are mostly using O I as an indicator (see JO} . 

The relative coverages of velocities where fl 
and /2 could be measured, weighted by their re- 
spective values of N(C II), are given for each sight- 
line in Column (4) of the table. These quantities 
vary by large factors from one case to the next. 
For the entire survey that spanned a total length 
of 180 kpc, the total N(C II) = 3.8 x 10 19 ciTr 2 , 
while that within our sampled velocity intervals 
represents N(C II) = 2.3 x 10 19 c m - 2 (61%). It 
is difficult to gauge the real fraction of the gas 
for which our measurements apply (i.e., CNM vs. 
CNM + WNM) because some of the WNM ma- 
terial can overlap in velocity the CNM that we 
sampled. Based on approximate interpolations of 
the velocity profiles of gas that is relatively free 
of C I, we estimate the fraction to be in the gen- 
eral vicinity of 15%, which means that on average 
our determinations could be systematically low by 
about —0.07 dex. 

Uncertainties in <7i ow and log(p/k)\ ow are prob- 
ably dominated by deviations of the real high- 
pressure conditions from those that apply to our 
adopted location for /l,/2 = (0.38,0.49), as we 
outlined in $5] One can estimate the magnitudes 
of such deviations by examining plausible alterna- 
tive geometrical constructions of the type depicted 
in Fig. 21 For example, if conditions in the high 
pressure gas are closer to T = 100 K, log(p/fc) = 
5.3 and #high ~ 0.1 (i.e., twice the general aver- 
age), an apparent value of log(p/fc)i ow = 3.5 may 
be 0.1 dex higher than the true value. The magni- 
tude of this effect scales in proportion to <?high, and 
it is diminished for higher values of log(p/k)i ow . 

In the light of our remarks about various forms 
of sampling bias in <j2](item 1), it should come as 
no surprise that values of log(//Jo) shown in Col- 
umn (7) of the table are all greater than zero. This 
is a consequence of the CNM being preferentially 
located in the vicinity of hot stars, rather than in 
random locations in the Galactic disk. This prefer- 
ence seems to overcome the effects of attenuation 
of starlight by dust. However, to some limited ex- 
tent our intensity outcomes could be elevated in 
a systematic fashion by the presence of unrelated 
WNM gas that is at the same velocity as the C I. 
This extra gas would mislead us into thinking the 
carbon atoms in the regions of interest are more 
ionized than in reality. 
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Table 3 

Observed and Calculated Quantities over all Velocities in the Sightlines 



Target 




Ccilc. log 


Percent C II 


Wei p 


hted Averages 1 ^ 


fvlcdian 


Star 


(cm' 2 ) 


NIC ll) b (cm _2 l 


Obscrvcd c 


9\ow 


lu B\y/ "'/low 




I'll 


(2) 


(1) 


(41 


(51 


(6) 


(7) 
V 1 ) 


CPD-59D2603 


>14.64 


17.87 


25.4 


0.92 


3.63 


0.31 


HD 108 


§14.98 


17.81 


64.4 


0.95 


3.82 


0.30 


HD 1383 


§14.80 


17.84 


44.5 


0.96 


3.49 


0.27 


HD 3827 


13J57±0.07 


17.22 


7.3 


0.92 


3.51 


0.16 


HD 15137 


14.62±0.01 


17.69 


74.3 


0.98 


3.61 


0.32 


HD 23478 


>14.75 


17.35 


46.4 


0.93 


3.62 


0.46 


HD 24190 


14.45±0.01 


17.49 


84.5 


0.97 


3.64 


0.62 


HD 24534 (X Per) . . 


>14.81 


17.49 f 


5.0 


0.95 


4.17 


1.16 


HD 27778 (62 Tau) . . 


>14.98 


17.48 f 


53.4 


0.90 


3.54 


0.49 


HD 32040 


13.26±0.03 


16.58 


18.1 


0.97 


3.07 


0.11 


HD 36408 


>14.17 


17.33 


54.5 


0.93 


3.82 


0.81 


HD 37021 (e 1 Ori B) 


13J34±0.05 


17.78 f ' 


67.5 


0.39 


3.83 


1.83 


HD 37061 (v Ori) . . . 


13.97±0.03 


17.90 f 


86.3 


0.59 


4.28 


1.80 


HD 37903 


14.24±0.05 


17.75 


48.8 


0.79 


4.61 


1.37 


HD 40893 


14.67±0.01 


17.79 


91.7 


0.98 


3.40 


0.37 


HD 43818 (11 Gem) . 


14.76±0.01 


17.89 


91.9 


0.98 


3.49 


0.35 


HD 44173 


13.66±0.05 


16.98 


4.1 


0.83 


3.45 


0.45 


HD 52266 


14.34±0.01 


17.62 


63.0 


0.98 


3.39 


0.48 


HD 69106 


14.30±0.01 


17.37 


81.1 


0.96 


3.55 


0.45 


HD 71634 


14.13±0.02 


17.18 


68.3 


0.89 


4.08 


0.88 


HD 72754 (FY Vcl) . 


gl4.40 


17.54 


37.0 


0.89 


3.86 


0.83 


HD 75309 


14?39±0.02 


17.51 


74.3 


0.97 


3.41 


0.46 


HD 79186 (GX Vcl) . 


>14.57 


17.76 


69.8 


0.96 


3.33 


0.40 


HD 88115 


14J33±0.05 


17.56 


29.4 


0.98 


3.55 


0.51 


HD 91824 


§14.45 


17.49 


63.6 


0.93 


3.62 


0.78 


HD 91983 


14J54±0.01 


17.55 


76.3 


0.96 


3.53 


0.51 


HD 93205 


14.56±0.02 


17.83 


51.5 


. . . g . . . g 


0.58 


HD 93222 


14.36±0.01 


17.82 


72.6 


0.95 


4.41 


0.82 


HD 93843 


14.15±0.05 


17.67 


12.8 


0.88 


4.13 


0.66 


HD 94454 . . 


>14.29 


17.50 


30.5 


0.93 


3.60 


0.68 


HD 94493 


14.26±0.02 


17.51 


59.1 


0.97 


3.53 


0.49 


HD 99857 


14.59±0.01 


17.71 


69.9 


0.94 


3.65 


0.55 


HD 99872 


>14.41 


17.45 


79.4 


0.94 


3.62 


0.86 


HD 102065 


>14.22 


17.40 


53.9 


0.96 


3.62 


0.61 


HD 103779 


14^23±0.03 


17.60 


44.8 


1.00 


3.29 


0.33 


HD 104705 (DF Cm) 


14.25±0.01 


17.58 


66.7 


0.98 


3.48 


0.53 


HD 106343(DL Cru) 


14.34±0.01 


17.62 


58.1 


0.96 


3.51 


0.49 


HD 108002 


>14.45 


17.58 


51.3 


0.96 


3.41 


0.34 


HD 108639 


14?31±0.01 


17.75 


76.3 


0.96 


3.49 


0.55 


HD 109399 


14.27±0.02 


17.60 


37.9 


0.97 


3.69 


0.59 


HD 111934 (BU Cru) 


>14.48 


17.70 


43.5 


0.99 


3.71 


0.63 


HD 112999 


14^23±0.02 


17.42 


82.7 


0.97 


3.61 


0.53 


HD 114886 


14.73±0.01 


17.75 


66.5 


0.95 


3.62 


0.37 


HD 115071 


14.69±0.01 


17.87 


74.3 


0.94 


3.64 


0.74 


HD 115455 


14.63±0.02 


17.85 


61.3 


0.97 


3.52 


0.58 


HD 116781 


14.28±0.02 


17.75 


32.2 


0.98 


3.42 


0.47 


HD 116852 


14.15±0.02 


17.45 


66.0 


0.95 


3.72 


0.80 


HD 120086 


13.20±0.08 


16.98 


3.6 


0.99 


3.72 


0.50 


HD 121968 


13.38±0.06 


16.87 


47.0 


0.95 


3.81 


1.27 
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Table 3— 


Continued 








Target 


logivrc i tot;i i) a 

o \ total / 


Calc. log 


Percent C II 


Weighted Averages' 1 


Median 


Star 


(cm -2 ) 


N(C II) b (cm 


Obscrved c 


Slow 


log(p/fc)iow 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


HD 122879 


14.42±0.01 


17.75 


76.3 


99 


3.59 


0.49 


HD 124314 


§14.66 


17.85 


66.0 


0.97 


3.54 


0.58 


HD 140037 


§14.03 


17.24 


4.8 


pa 1.0 


5.19 


0.95 


HD 142315 


13?78±0.06 


17.33 


7.3 


0.93 


3.38 


0.50 


HD 142763 


13.38±0.06 


17.03 


24.9 


0.93 


3.47 


0.70 


HD 144965 


>14.28 


17.52 


9.4 


0.82 


3.80 


0.76 


HD 147683 


§14.80 


17.64 


39.7 


0.89 


3.89 


0.57 


HD 147888 (p Oph) 


>14.51 


17.80 f 


26.9 


0.68 


3.98 


1.20 


HD 148594 


§14.12 


17.54 


39.3 


n 7Q 
u. / y 


A M 
'■±.01. 


1.53 


HD 148937 


§14-87 


17.96 


90.6 


u.yu 


3 RA 


0.67 


HD 152590 


14^61±0.01 


17.76 f 


77.0 


U.O 1 


1 71 
O. ( 1 


0.77 


HD 156110 


13.88±0.04 


17.08 


45.6 


0.99 


3.83 


0.57 


HD 157857 


14.61±0.01 


17.77 


87.2 


u.yu 


o.ou 


0.46 


HD 165246 


§14.33 


17.73 


59.4 


0.94 


3.58 


0.79 


HD 175360 


13?98±0.01 


17.29 


76.6 


1 no 


•3 1Q 
O. lO 


0.39 


HD 177989 


§14.66 


17.44 


72.5 


u.yu 


o.oy 


0.35 


HD 185418 


>14.72 


17.71 


66.8 


0.97 


3.41 


0.23 


HD 192639 


§14.74 


17.85 


80.0 


n Q7 
u.y i 


o.uo 


0.52 


HD 195965 


§14.48 


17.42 


60.2 


u.yu 


o.ou 


0.32 


HD 198478 (55 Cyg) 


§14.84 


17.80 


84.5 


n Q9 
u.yz 


o.uo 


0.43 


HD 198781 


§14.56 


17.57 


32.0 


Q9 
u.yz 


o.^y 


0.37 


HD 201345 


13^97±0.02 


17.42 


63.4 


0.97 


3.50 


0.32 


HD 202347 


14.61±0.01 


17.33 


79.5 


0.95 


3.76 


0.20 


HD 203374 


§14.98 


17.68 


81.1 


0.94 


3.63 


0.31 


HD 203532 


>14.65 


17.40 


4.2 


0.85 


4.40 


1.26 


HD 206267 


§15.29 


17.85 


72.4 


0.93 


3.64 


0.30 


HD 206773 


§14.70 


17.55 


71.8 


0.96 


3.55 


0.19 


HD 207198 


§15.24 


17.81 f 


87.6 


0.94 


3.63 


0.20 


HD 208440 


§14.84 


17.68 


78.8 


0.94 


3.66 


0.35 


HD 208947 


§14.63 


17.36 


53.2 


0.97 


3.60 


0.33 


HD 209339 


14?76±0.01 


17.61 


87.1 


0.96 


3.69 


0.34 


HD 210809 


§14.70 


17.69 


36.8 


0.94 


3.66 


0.29 


HD 210839 (A Cep) . 


§14.95 


17.79 


63.8 


0.86 


4.16 


0.47 


HD 212791 


14?28±0.03 


17.44 


29.9 


0.97 


3.54 


0.38 


HD 218915 


14.56±0.01 


17.59 


68.9 


0.97 


3.58 


0.33 


HD 219188 


13.92±0.04 


17.08 


77.6 


0.99 


2.97 


0.01 


HD 220057 


>14.71 


17.46 


25.6 


0.94 


3.51 


0.35 


HD 224151 


§14.62 


17.80 


56.1 


0.97 


3.80 


0.16 


HDE 232522 


§14.60 


17.66 


58.5 


0.97 


3.53 


0.35 


HDE 303308 


14?69±0.00 


17.83 


72.9 


... g 


... g 


0.59 



integrated over all velocities where C I absorption seems to be visible (not just over the restricted regions 
where the lines are strong enough to yield good measurements of /l and /2). Sometimes there was evidence that 
unresolved saturations were evident at certain velocities, as indicated by a test that is discussed in ^A.2.41 When 
this occurred over very limited portions of the profile, we indicate a mild inequality with "§." When a substantial 
portion of the profile exhibited such behavior, we indicate a more severe inequality by ">." When errors in the 
column densities are given, they indicate only the quantifiable errors arising from noise or uncertainties in the 
continuum levels. These errors indicate the relative quality of the measurements, but they are not fully realistic 
because they do not take into account uncertainties in our adopted /-values or occasional flaws in the MAMA 
detector used by STIS. 

b The computed amount of C II at all velocities based on the absorption profiles of O I or S II; see i|4,ll for 
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details. These amounts compare favorably with the observed amounts for a few stars in note /. 

c The relative amount of C II, as represented by its proxy O I (and sometimes S II), within the velocity interval 
where determinations of fl and /2 were good enough to be considered for pressure measurements, compared to 
the amount seen at all velocities, as shown in the previous column. 

d Calculated according to the following: J2[9lowN(C I to tal)]/ N(C Itotal) and 

log£[(p/fc) low AT(C II)]/£iV(C II). 

e Our estimate for the local density of UV radiation from starlight that is more energetic than the ionization 
potential of neutral carbon, compared to an adopted standard Iq based on a level specified by Mathis et al. (1983) 
for the average intensity of ultraviolet starlight in the Galactic plane at a Galactocentric distance of 10 kpc. This 
estimate is based on our evaluation of the ionization equilibrium of C, as expressed in Eq. [l] 

f Compare with actual measurements of logiV(C II) using the intersystem C II] line at 2325 A: From 
Sofia et al. (1998) HD 24534: 17.51 (+0.11, -0.16). From Sofia et al. (2004) HD2778: < 17.34; 
HD 37021: 17.82 (+0.12, -0.18); HD 37061: 18.13 (+0.04, -0.06); HD 147888: 18.00 (+0.07, -0.09); HD 152590: 
18.21 (+0.08, -0.10); HD 207198: 17.98 (+0.11, -0.14). However, recent measurements of the damping wings for 
the allowed transition at 1334.53 A reported by Sofia et al. (2011) indicate that these column densities may be 
too large by a factor of about 2. 

g Gas within a component at large negative velocities has conditions very near the high pressure reference mark. 
Hence the projection onto the low pressure arc is meaningless. 
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Of particular interest are the characteristic sizes 
of the regions containing the C I that we are able to 
study. Within any velocity bin, we can determine 
a value for the local density of gas particles, com- 
posed of atomic hydrogen, helium atoms, and hy- 
drogen molecules. Once again, if we assume that 
/(H 2 ) = 0.6 (see JO]) and He/H = 0.09, it follows 
that the local density of hydrogen nuclei is given 
byn(H) = p/(0.79kT). The longitudinal thickness 
occupied by the gas is then equal to the column 
density of these nuclei, -/V(H), divided by n(H). We 
obtain -/V(H) by multiplying the amount of carbon, 
measured by the methods outlined in §4.11 by the 
general expectation for (H/C) = 5040 in the ISM. 
A sum over velocity of all of the length segments 
gives the overall thickness of the C I-bearing re- 
gion(s) in any particular line of sight. 

Figure O shows the distribution of region thick- 
nesses for all of the lines of sight in our survey, 
which generally have dimensions of less than 20 pc. 
The occupation fractions are quite small, generally 
less than about 2%. 

For the benefit of future investigations that 
may require more detailed information about in- 
dividual sight lines, Table @] provides a machine- 
readable summary that lists for each velocity chan- 
nel the measured values of fl, /2, and N(C hotai), 
along with the calculated values of N(C II), g\aw, 
log(>/fc) low , and log(I/I ). 
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Fig. 5. — Histograms that show (left) the total 
thicknesses of the regions and (right) their relative 
occupation fractions in the sightlines that we are 
able to measure in the survey. 



8. Behavior with Velocity 

The radial velocities that we measure in the C I 
profiles arise from various kinematical phenomena, 
such as differential velocities caused by rotation 
or density waves in the Galaxy, coherent motions 
caused by discrete dynamical events such as super- 
nova explosions, mass loss from stars, the collision 
of infalling halo gas with material in the Galactic 
plane, and random motions arising from turbu- 
lence. With the exception of differential Galactic 
rotation, all of these effects can transform some of 
their energy into an increase of the thermal pres- 
sures. In their limited sample of only 21 stars, 
JT01 found elevated pressures in gases whose ve- 
locities deviated away from the range that was ex- 
pected for differential Galactic rotation. We now 
revisit this issue for our present, much larger sam- 
ple of sightlines to further substantiate the evi- 
dence for a coupling between the thermal pressures 
and unusual dynamical properties of the gas. 

Figure |5] shows the measurements of /l and f2 
at velocities that are either above or below the re- 
spective line-of-sight velocity ranges permitted by 
differential Galactic rotation, assuming that the 
rotation curve is flat at 254 km s^ 1 and the dis- 
tance to the Galactic center is 8.4 kpc (Reid et 
al. 2009). An extra margin of 5km s _1 is added 
to the exclusion zone for permitted velocities, so 
that we are more certain of showing material that 
genuinely disturbed in some manner. When we 
compare the results of Fig. [6] to those shown in 
Fig. [5J it is clear that gases at high velocity do 
not have nearly as strong a central concentra- 
tion near fl w 0.2 and /2 w 0.07. The "center 
of mass" (fl, f2) locations for all of the points 
shown in Fig. [6] are (0.265, 0.163) for negative ve- 
locities and (0.228, 0.078) for positive velocities. 
By comparison, for all measurements shown in 
Fig. [21 we found the balance point to be at (0.209, 
0.068). These differences are principally caused by 
a greater prominence of a more highly dispersed 
population of points in fl and /2, and they should 
come as no surprise since they demonstrate the ex- 
pected coupling of the dynamics of the gas to the 
observed enhancements in the thermal pressures. 
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Fig. 6. — Presentations similar to Fig. [2j except that the measurements include velocities only below the 
minimum value permitted by differential Galactic rotation, but with an extra margin of 5km s — , i.e., 
v < min(u gr ) — 5 km s _1 (left-hand panel) or more than 5 km s _1 above the maximum permitted velocities, 
i.e., v > max(wg r ) + 5km s -1 (right-hand panel). The dot diameters in these diagrams are twice as large as 
those in Fig. [5] for a given column density of C Itotai- 
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Table 4 



Observed and Calculated Quantities in Specific Velocity Channels 


Target 


Velocity 


/I 


/I 


/2 


/2 


N(C I total ) 


N(C II) 


9low 


log(p/fc) low 


log(/// ) 


Star 


(km s -1 ) 




Error 




Error 


(cm -2 ) 


(cm -2 ) 








CPD-59D2603 


4.0 


0.213 


0.016 


0.068 


0.026 


2.13c+013 


1.34c+016 


0.930 


3.63 


0.37 


CPD-59D2603 


4.5 


0.201 


0.013 


0.065 


0.020 


2.81c+013 


1.44e+016 


0.933 


3.60 


0.28 


CPD-59D2603 


5.0 


0.200 


0.011 


0.065 


0.017 


3.33e+013 


1.56e+016 


0.934 


3.60 


0.25 


CPD-59D2603 


5.5 


0.207 


0.010 


0.065 


0.015 


3.83e+013 


1.76e+016 


0.936 


3.62 


0.26 


CPD-59D2603 


6.0 


0.218 


0.009 


0.064 


0.013 


4.46e+013 


2.02e+016 


0.945 


3.67 


0.29 


CPD-59D2603 


6.5 


0.224 


0.008 


0.065 


0.011 


5.30e+013 


2.46e+016 


0.944 


3.68 


0.31 


CPD-59D2603 


7.0 


0.227 


0.008 


0.066 


0.009 


6.26e+013 


3.03e+016 


0.943 


3.69 


0.33 


CPD-59D2603 


7.5 


0.230 


0.008 


0.070 


0.008 


7.17c+013 


3.80e+016 


0.938 


3.70 


0.37 


CPD-59D2603 


8.0 


0.233 


0.008 


0.074 


0.007 


7.88e+013 


4.63e+016 


0.931 


3.70 


0.41 


CPD-59D2603 


8.5 


0.222 


0.008 


0.079 


0.007 


8.10c+013 


5.13c+016 


0.910 


3.65 


0.40 


CPD-59D2603 


9.0 


0.200 


0.007 


0.081 


0.007 


7.65c+013 


5.18e+016 


0.894 


3.56 


0.35 


CPD-59D2603 


9.5 


0.181 


0.007 


0.081 


0.009 


6.53e+013 


4.85e+016 


0.885 


3.47 


0.32 


HD102065 . . . 


8.5 


0.188 


0.019 


0.071 


0.016 


2.42c+012 


5.98e+014 


0.908 


3.51 


0.09 


HD102065 . . . 


9.0 


0.241 


0.012 


0.072 


0.010 


3.87e+012 


6.54e+014 


0.934 


3.71 


0.14 


HD102065 . . . 


9.5 


0.275 


0.008 


0.067 


0.006 


5.93e+012 


7.52c+014 


0.966 


3.83 


0.14 


HD102065 . . . 


11.0 


0.223 


0.004 


0.050 


0.003 


1.64c+013 


1.56e+016 


0.973 


3.62 


0.61 


HD102065 . . . 


11.5 


0.203 


0.005 


0.048 


0.002 


1.92e+013 


1.99e+016 


0.967 


3.55 


0.59 


HD102065 . . . 


12.0 


0.187 


0.005 


0.046 


0.002 


2.14c+013 


2.40e+016 


0.964 


3.49 


0.57 


HD102065 . . . 


12.5 


0.176 


0.005 


0.043 


0.002 


2.38e+013 


2.83e+016 


0.965 


3.44 


0.54 


HD102065 . . . 


13.0 


0.177 


0.004 


0.043 


0.002 


2.70e+013 


3.13e+016 


0.965 


3.45 


0.55 


HD102065 . . . 


13.5 


0.187 


0.004 


0.046 


0.002 


3.06e+013 


3.67c+016 


0.967 


3.49 


0.59 


HD102065 . . . 


14.0 


0.204 


0.004 


0.051 


0.002 


3.32c+013 


4.35e+016 


0.964 


3.54 


0.66 


HD102065 . . . 


16.0 


0.318 


0.005 


0.099 


0.003 


1.47c+013 


2.83c+016 


0.925 


3.85 


1.06 


HD102065 . . . 


16.5 


0.315 


0.005 


0.101 


0.004 


1.00e+013 


1.87e+016 


0.918 


3.84 


1.03 


HD102065 . . . 


17.0 


0.298 


0.007 


0.091 


0.006 


6.78e+012 


1.05c+016 


0.931 


3.81 


0.95 


HD102065 . . . 


17.5 


0.281 


0.010 


0.074 


0.008 


4.54e+012 


3.97e+015 


0.954 


3.80 


0.74 


HD102065 . . . 


18.0 


0.270 


0.014 


0.059 


0.012 


3.07e+012 


2.74c+015 


0.983 


3.78 


0.72 


HD102065 . . . 


18.5 


0.273 


0.021 


0.053 


0.018 


2.11c+012 


2.55e+015 


0.998 


3.77 


0.80 


HD102065 . . . 


19.0 


0.283 


0.030 


0.055 


0.026 


1.45c+012 


2.15c+015 


1.000 


3.79 


0.89 



(This tabic is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 
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Models of pressurized clouds behind weak 
shocks in the ISM computed by Bergin et al. 
(2004) reveal that the outcomes for (fl, /2) are 
centered on values of approximately (0.36,0.18), 
(0.40, 0.25) and (0.37, 0.35) for the post-shock con- 
densations behind shocks with velocities of 10, 
20 and 50 km s~ , respectively (see their Fig. 8; 
in these cases the resultant ram pressures were 
1.4 x 10 4 , 5.8 x 10 4 and 3.6 x 10 5 citT 3 K for a 
preshock density of 1cm -3 ). These results for fl 
and /2 are well removed from the densest cluster- 
ing of measurements shown in Fig. [2j but they do 
seem consistent with the more sparse population 
of points having f2 > 0.15, which is more strongly 
emphasized in the unusual velocity ranges repre- 
sented by the two panels of Fig|6l 

Figure [5] shows clearly that a moderate num- 
ber of the measurements in the positive-velocity 
regime exhibit higher pressures than usual, but 
not to the great extremes revealed by the nega- 
tive velocity gas. We offer a simple interpretation 
for why this happens. We propose that a signifi- 
cant fraction of the high pressure material arises 
from stellar mass-loss outflows that eventually col- 
lide with the ambient medium, creating dense, ex- 
panding shells that are at high pressures (Castor 
et al. 1975; Weaver et al. 1977). Another possibil- 
ity is that small clouds surrounding the stars are 
pressurized and accelerated by either the momen- 
tum transfer arising from photoevaporation (Oort 
& Spitzer 1955; Kahn 1969; Bertoldi 1989) or the 
momentary surge in pressure of a newly developed 
H II region. These phenomena associated with our 
target stars should be visible to us if they contain 
C I. The shells should be intercepted by our sight 
lines regardless of whether they (or possibly small 
clouds inside them) are large and at moderately 
high pressures or very small and at much higher 
pressures. The foreground portions of such shells 
are responsible for the negative velocity gas that 
we can view in the star's spectrum. For positive 
velocity gas the situation is different. Here, we 
rely entirely on the random chance of seeing ei- 
ther one of the large-scale events (of non stellar 
origin) mentioned at the beginning of this section, 
or else the rear portion of some region or shell 
that is created by some foreground star or stellar 
association that is unrelated to the star that we 
are viewing. This being the case, there may be 
a vanishingly small chance that we will intercept 



a highly pressurized shell with a small diameter, 
but the chances increase for larger shells that have 
lower pressures at their boundaries. This obser- 
vational bias against small, high pressure events 
at positive velocities could explain why we see no 
points at \og(p/k) > 5 in the right-hand panel of 
the figure. 

In the next section (and in §10.5. 4j) . we will re- 
inforce the picture that high pressures arise from 
the increased dynamical activity near bright stars. 
We will show evidence that there is a strong corre- 
lation between pressures and the local intensities 
of ionizing radiation. 

9. Interpretations of the Results 
9.1. Basic Distribution Functions 

After evaluating the conditions within each ve- 
locity interval for all of the lines of sight, we are 
in a position to look at the composite outcome 
of all of the results of the dominant low-pressure 
component. All of the presentations in this sec- 
tion will show distributions expressed in terms of 
the amount of hydrogen in a given condition. In 
order to do so, we must convert our original mea- 
surement weights based on N(C Itotal) into ones 
that account for the equivalent column densities 
of C II that we derived from our determinations of 
O I (and on rare occasions S II) at identical veloc- 
ities, as discussed in i j4.ll Once again, we convert 
from N(C II) to iV(H) by multiplying the amount 
of carbon by (H/C) = 5040 in the ISM. As we in- 
dicated earlier ( §4.4|) . WNM material at the same 
velocity as the CNM will tend to inflate somewhat 
the derived value of iV(H) associated with the C I 
that is used for determining the pressure. 

The histogram distributions shown in Fig. [7] re- 
veal that the pressure distribution function is not 
strongly influenced by the temperature of the gas, 
as deduced from the measurements of T i of H 2 . 
Nevertheless, the evidence that we have suggests 
an inverse correlation of pressures with tempera- 
ture, although the scatter in this relationship is 
large. For the points shown in the inset of the 
figure, the Spearman rank order correlation coeffi- 
cient is —0.29. This determination is significantly 
different from a zero correlation for the popula- 
tion at the 97.5% confidence level for 58 pairs of 
measurements. The dispersion of the results is so 
large that it is difficult to assign a value for the 
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Fig. 7. — The distribution of thermal pressures, normalized to the estimated amount of hydrogen present, for 
three different temperature conditions, as indicated by the J = to 1 rotation temperatures of H2: (1) all of 
the gas (tallest profile), (2) gas for which Tqi > 85 K (middle profile) and (3) gas for which T < 75 K (shortest 
profile). The median temperature for all cases is 77 K, i.e., a value that is between the two limits. Sight lines 
where Tqi measurements do not exist are included in condition (1) but excluded from conditions (2) and 
(3). A best-fit lognormal distribution for condition (f) is shown by the solid curve, and it is represented by 
Eq. [3] The inset shows a scatter plot of the C I weighted average \og(p/k) given in Column (6) of Tabled 
vs. T i (if known), as listed in Column (7) of Table [5J 



28 



o 



o 




\og(p/k) [cm" J K] 



Fig. 8. — A gray-scale representation of the loga- 
rithm of the amount of H I gas that we found as a 
function of log(J/Jo), i- e -> the logarithm of the en- 
hancement of the starlight density above average, 
against the thermal pressure, expressed in terms of 
log{p/k). The dashed line shows the cutoff equal 
to y/lO times the average field that we established 
for defining the low intensity distribution shown 
in Fig. H 

apparent polytropic index of the gas, but the sign 
of the trend is consistent with a slope of less than 
one for the CNM thermal equilibrium track near 
the minimum pressure shown in Fig. [TJ 

The central portion of the distribution of ther- 
mal pressures (for all Toi) follows closely a lognor- 
mal distribution given by 



dN(R)/dlog(p/k) 



2.30 x 10 23 exp 



(log(p/fc)-3.58) 2 
2(0.175)^ 



(3) 



This lognormal relationship is shown by the 
smooth curve in Fig. [7] (and will be shown 
again later in a log-log representation by the 
smooth gray curve in Fig. [S]). Outside the range 
3.2 < log(p/fc) < 4.0 it understates the observed 
amount of material in the wings of the profile (this 
is not evident in Fig. [71 but is clearly shown later 
in Fig.i]). 

Figure [S] shows that the outcomes for the 
starlight densities and the thermal pressures are 



not independent of each other. In regions that are 
close to stars that emit UV radiation (//Jo ~ 10), 
we find that with a few exceptions the average 
pressures increase to values in the general vicinity 
of log(p/fe) ~ 4. This enhancement supports the 
viewpoint that turbulent energies are greater in 
the general vicinity of young stars, a phenomenon 
that may be related to changes in the morphology 
of H I near stellar associations that were found by 
Robitaille et al. (2010). 

In order to obtain a representation for the pres- 
sures in the general ISM somewhat removed from 
the bright stars, we will limit further study of 
the distribution to only those cases where 1/ Jo < 
10 5 , a limit that is depicted by the dashed line in 
Fig. [8] We consider that any gas elements that are 
above that line represent localized regions that are 
exceptionally close to sources of mechanical energy 
and are thus not representative of the general, dif- 
fuse ISM. 

Figure [9] shows in a log-log format the dis- 
tribution of thermal pressures for J/Jq < 10 5 
(black histogram) compared with the distribu- 
tion for all intensities (gray histogram). In terms 
of log(p/fc) [and using a linear representation of 
dN(H)/d\og(p/k)], the distribution for the low- 
intensity results has a mean of 3.47, a standard 
deviation of 0.253, a skewness of —1.8, and a kur- 
tosij"! of 6.2. The influence of the excess of low 
pressure outcomes, as evidenced by the negative 
skewness, causes the standard deviation listed here 
to be larger than the value 0.175 given in Eq. [3l 
which would apply to just the central portion of 
the profile. 

For the convenience of those who wish to re- 
produce a reasonably good representation of the 
low-intensity data in analytical form, we supply 
an empirical polynomial fit, 

dN(U)/dlog(p/k) = 1.16 x 10 23 exp(-0.0192z 



-0.00387z 2 
-6.24 x 10" V - 6.77 x 10 



1-5,3 



-9,5 



z ) cm 



(4) 



where the dimensionless quantity z — (p/k)z — 60 
(for p/k expressed in terms of cm~ 3 K). This em- 



10 Our definition of kurtosis includes a subtraction of 3 from 
the fourth moment divided by <r 4 , thus making the kurtosis 
of a Gaussian distribution equal to zero. Sometimes in the 
literature, e.g. Federrath et al. (2010), this "—3" term is 
omitted. 
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Fig. 9. — Log-log presentations of the distribution of thermal pressures for two cases: all of the gas sampled 
by C I is shown by the gray histogram, while a subset of the material for which I/Iq < 10° 5 is shown by 
the black histogram. This intensity cutoff limits the sample to all of the outcomes that appear below the 
dashed line in Fig. [8j The thin curves show how well the analytical expressions given in Eqs. [3J (gray) and 
|4] (black) fit the results. The vertical dot-dash line labeled CNM \og(p/k) mm corresponds to the minimum 
pressure that is allowed for a static CNM, as shown by a similar line with the same designation in Fig. [T] 



pirical fit is shown by the thin, black curve in 
Fig. |U If this distribution function is converted 
into a linear representation, i.e., AT(H) as a func- 
tion of p/k, we find that for p/k < 5500 cm~ 3 K it 
docs not differ appreciably from a Gaussian func- 
tion with mean value of p/k = 3700 cm~ 3 K and a 
standard deviation of 1200cm~ 3 K. The distribu- 
tion is somewhat higher than this Gaussian func- 
tion above 5500 cm -3 K. 

Our mean value quoted above is 0.22 dex higher 
than the value 2240cm _3 K that we listed earlier 
(JT01). There are three independent reasons that 
can account for nearly all of this difference. First, 



we used revised rates for the collisional excitation 
and radiative decay of the upper fine-structure lev- 
els of C I, as discussed in E1A.3I This accounts 
for an elevation of typical determinations of pres- 
sures of about 0.05 dex. Second, our new esti- 
mate for the strength of the optical pumping of 
the levels has been reduced ( §A.4p . with the re- 
sult that a typical pressure measurement should 
be raised by another 0.05 dex. Third, our earlier 
specification for the mean value of p/k was for a 
temperature (40 K) that gave the lowest inferred 
pressure for a given level of C I excitation, while 
the present result uses either actual temperatures 



30 



measured from the H2 rotational excitations or a 
median value of 80 K if an explicit measurement 
for a sight line is not available. Under most cir- 
cumstances, the inferred pressure at 80 K is about 
0.1 dex higher than that for 40 K. Taken together, 
these three effects can account for an elevation of 
our new pressures over the old ones by 0.2 dex. 

As a final point, we add a cautionary note that 
the errors in our determinations for log(p/fc) be- 
come much larger than usual when their values fall 
below 3.0. The reason for this is that the changes 
in fl, the major discriminant for pressures, be- 
come very small at low pressures, as shown by the 
shrinkage in the spacing between the 0.1 dex mark- 
ers in Figs. [5] and |U 

9.2. Volume- Weighted Distributions 

Up to now, the distribution functions that we 
have shown (Figs.[7H§]) have been weighted in pro- 
portion to our calculated hydrogen column densi- 
ties, which is equivalent to a sampling by mass. In 
many cases, investigators showing results of com- 
puter simulations of ISM turbulence express their 
outcomes according to the counts of volume cells 
having different pressures. In order to make a 
conversion from a mass-weighted distribution to 
a volume- weighted one, we must make a simplify- 
ing assumption that we are viewing an ensemble 
of gases that has internal random pressure fluctua- 
tions that change with time, but that is otherwise 
approximately uniform in nature and that can be 
characterized as having an equation of state with 
a single value for the polytropic index 7 (equal to 
the slope of logp vs. logn or the ratio of specific 
heats Cp/c v ). In this situation, the changes in pres- 
sure cause the volumes of mass parcels to change 
in proportion to p _1//7 , and this factor must be ap- 
plied to the mass- weighted distribution function to 
obtain the volume-weighted one. 

The smooth curves in Fig. [TU] show how the 
mass- weighted distribution in \og(p/k) for the low 
starlight intensities would appear after being con- 
verted to volume weighted ones for three differ- 
ent assumed values of 7. These three examples 
illustrate the behavior of the gas under the con- 
ditions (1) 7 = 0.7, which is a good approxima- 
tion of the slope of thermal equilibrium curve for 
the CNM shown in Fig. [TJ (2) the relationship for 
7 = 1.0 that corresponds to an isothermal gas, and 
(3) a condition 7 = 5/3 that indicates that the 



gas is undergoing adiabatic compressions and ex- 
pansions (and assuming that the gas has a purely 
atomic composition). The divergent behavior of 
the curves at the far left portion of this diagram 
probably arises from either deviations caused by 
small number statistics for the samples at the ex- 
tremely low pressures or the fact that the errors in 
log(p/k) become larger than usual at the low pres- 
sure extreme. It is important to emphasize that 
in a turbulent cascade the notion that the gas has 
a single polytropic index on all length scales is 
an oversimplification; we will explore this issue in 
more detail in §10.31 

9.3. Pileup in Velocity Bins 

As we discussed in Ej3j any outcomes for fl and 
/2 at a particular velocity may represent a com- 
posite result for two or more regions that are seen 
in projection along the line of sight. In 35]we ex- 
plained how we separated contributions from small 
amounts of gas at extraordinarily large pressures, 
well away from the dominant regime of low pres- 
sures. However, we have yet to address the possi- 
bility that two or more regions at somewhat dif- 
ferent pressures along the lower, nearly straight 
portion of the fl — /2 equilibrium curve can cre- 
ate an apparent outcome that represents a proper 
C I-weighted mean, but without revealing the true 
dispersion of pressures from the contributors. If 
such superpositions are taking place frequently, 
they will tend to decrease the width our observed 
overall distribution shown in Figs. [THU 

One way to gain an insight on this possibil- 
ity is to examine how deviations from the mean 
log(p/fc) scale with the corresponding amounts of 
C I. If we imagine a simple picture where all of the 
C I exists within independent parcels, each with 
some single, representative value N (C Itotai), we 
would expect to find that the dispersion of any 
collection of measurements having some multiple 
n times -/V (C I to tai) would show us a standard de- 
viation equal to otruc/v^, where <7t rU c is the real 
dispersion in log(p/fc) for the individual packets 
that are seen in projection. 

The upper panel of Figure [11] shows the appar- 
ent outcomes for values of \og(p/k) as a function 
of iV(C Itotai) -0 ' 5 ; h is clear that as the column 
densities decrease (i.e., moving toward the right 
of the plot), the vertical dispersions increase. For 
data segregated within successive bins having a 
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Fig. 10. — (three smooth curves:) The results of a conversion of the mass-weighted distribution curve for 
starlight intensity levels I/Iq < 10° 5 (the black, smooth curve shown in Fig. |9]) to volume weighted ones for 
three assumed values for 7, corresponding to cases where the gas behaves as if it were in thermal equilibrium 
(7 = 0.70), isothermal (7 = 1.0) and adiabatic (7 = 1.67). (Histogram- style trace:) The volume-weighted 
distribution for all intensity levels for log(p/fc) > 3.5, assuming 7 = 1.67. This distribution is relevant 
to a discussion that is presented in M10.4I about the possible creation of higher than normal pressures by 
expanding supernova remnants. In all four cases, the curves are normalized such that their integrals over all 
\og(p/k) equal 1.0. 



width of V2 x 10~ 7 cm, the lower panel indicates 
that the rms dispersion indeed seems to scale in 
direct proportion to iV(C Itotai)~ ' 5 , but only up 
to about N(C Itotai) -0 5 = 5V2 x lO^cm (indi- 
cating that iV (C Itotai) ~ 2 x 10 12 cm- 2 ). Thus, 
on the one hand, one could imagine that <7tmc 
could be as large as around 0.5, instead of our 
overall measured value of 0.253. On the other 
hand, the proposed model for superpositions may 
be only a product of our imagination: perhaps 
coherent regions with larger values of N(C Itotai) 
have a real tendency be less easily perturbed by 
various external forces that cause pressure devia- 
tions away from some mean value. In essence, the 



trend shown in Fig. [TT] may reflect a real physi- 
cal effect rather than a trend caused by random 
superpositions of unrelated, small gas clouds. 

In principle, a trivial explanation for the ef- 
fect shown in Fig. FPU might be that as N(C Itotai) 
decreases the measurement errors in log(p/k) in- 
crease. However, as we explain later in EjA.2.31 
the la errors in /l and /2 should be equal to 0.03 
or less for the measurements to be accepted. At 
normal pressures this amount of error is equiva- 
lent to a change in \og(p/k) equal to only 0.1 dex, 
far smaller than the observed dispersion that is 
shown for low column density cases in Fig. [IT] One 
reason that we are able to maintain small errors 
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Fig. 11. — Top panel: Individual measurements of log(p/fc) in the low pressure regime as a function of 
the inverse square-root of the column density of C Itotai for all velocity channels of width 0.5 km s _1 that 
had log//io < lO 5 . The mean values of log(p/fc) are listed near the top of the panel for successive 
intervals centered on integral values of [N(C Itotai)/5 x 10 13 cm~ 2 ]~ ' 5 . Lower panel: The standard deviations 
in log(p/k) for measurements within the intervals, showing an almost linear progression up to a value 
er[log(p/fc)] f=a 0.5, at which point the column density of independent packets of gas have characteristic values 
N (C hotai) = 2 x 10 12 cm~ 2 . The "n =" designations show the number of points that were used to evaluate 
cr[\og(p/k)] in each bin. 



for lower column densities is that our system of 
weighting the measurements causes a shift of em- 
phasis from weak atomic transitions to stronger 
ones as iV(C Itotai) decreases. 

10. Discussion 

10.1. Distribution Width and Shape 

10.1.1. Overall Shape 

The highly symmetrical appearance of our dis- 
tribution for all of the material that we sampled in 
the regime of ordinary pressures (which we identi- 
fied as the "low pressure component" in Ej5j) is an 
illusion that arises from the projection of the irreg- 
ularly shaped distribution depicted in Fig. [8] onto 
the a;-axis that represents log(p//c). The distribu- 



tion function reverts to one with a strong negative 
skewness when we limit our consideration to con- 
ditions where I/Iq < 10 5 (below the dashed line 
in Fig. [8]) . This behavior is inconsistent with tur- 
bulence in an isothermal gas, which should show 
a pure lognormal density (and pressure) volume- 
weighted distribution function (Vazquez Semadeni 
1994; Nordlund & Padoan 1999; Kritsuk et al. 
2007). 

10.1.2. Deviations to Low Pressures 

A substantial fraction of the material (29%) - 
that which is depicted to the left of the verti- 
cal dash-dot line in Fig. [9] - is detected at pres- 
sures below those permissible for a static CNM, 
(p/k) m i n = 1960cm~ 3 K, as defined by the "stan- 
dard model" for the thermal equilibrium curve 
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presented by Wolfire et al. (2003) that we show 
in Fig. HHH From this we conclude that either 
their curve does not apply to the media we are 
viewing or else that rarefactions caused by tur- 
bulence create momentary excursions below the 
curve. The recovery toward normal pressures for 
regions that reach anomalously low pressures in 
some cases might be inhibited by temporary, lo- 
cally high values of magnetic pressure (Mac Low 
et al. 2005). 

A valid question to pose is whether or not we 
could understand the existence of low pressures by 
large changes in some of the parameters that in- 
fluence the value of (p/k) m i n within some localized 
regions. Wolfire et al. (2003) expressed a sim- 
ple equation [their Eq. (33)] that gives some guid- 
ance on this possibility. We restate their equation 
terms of our variables by substituting 0.674(//7o) 
for their normalized ISM intensity G' at a Galac- 
tocentric distance of 8.5 kpc. The reason for this 
substitution is that we have adopted a more re- 
cent measure of a standard intensity Iq (Mathis et 
al. 1983) that is lower than the one they chose to 
use [taken from Draine (1978)]. Also, we set their 
parameter for the total ionization rate (multiplied 
by 10 16 s _1 ) Ct = 2.0, since we have adopted a cos- 
mic ray ionization rate £cr = 2 x 10~ 16 s _1 (see 
^4.4p and assumed that the x-ray and EUV ion- 
ization rates are very small in comparison. Our 
restatement of their equation takes the form, 

(p/*)^ = 5730(^I/Io) 1 + 2 08( g J/Jo)0 , 365 , 

(5) 

where Z' d is equal to the normalized ratio of in- 
terstellar dust grains to polycyclic aromatic hy- 
drocarbons (PAH), and Z' g is the normalized gas 
phase abundance of heavy elements that are re- 
sponsible for radiative cooling of the gas (chiefly 
C and O). The two quantities Z' d and Z' g are gen- 
erally assumed to be equal to 1.0 for conditions 
in our part of the Galaxy. If log(//io) = —0.35, 
we find from Eq. [5] that (p/k) nlin = 1000cm~ 3 K. 
However, the distribution of outcomes shown in 
Fig. [S] indicate most of our pressure measurements 
apply to regions with log(J//o) > —0.35. Another 
way to reduce (p/k) m i n to 1000cm~ 3 K would be 



With a parametric formulation discussed in the next para- 
graph using our value for I/Iq = 1.0 and Ccr = 2 X 
10 — 16 s — 1 , (p/fc)min decreases very slightly to 1860 cm - 3 K. 



to have I/Iq = 1.0 but with either a ratio of the 
dust grain to the PAH concentration Z' d as low as 
0.45 times the normally assumed value or a reduc- 
tion of Z' g to 0.54 times the normal amount. Even 
with the possible deviations discussed here that 
would make (p/fc)min reach as low as 1000 cm _3 K, 
we still have measurable amounts of gas below the 
pressure threshold for a stable CNM. 

10.1.3. Comparisons with Expectations of the Ef- 
fects of Turbulence 

The magnitude and skewness of the fluctua- 
tions in thermal pressure give an indication of the 
strength and character of the turbulence in the 
ISM (Padoan et al. 1997b). For instance, the 
one-dimensional simulations of Passot & Vazquez- 
Semadeni (1998) illustrated how 7 changes the 
sign of the skewness of the distribution: 7 < 1 
makes the distribution shallower on the high pres- 
sure side of the peak and steeper on the low pres- 
sure side, while the opposite is true for 7 > 1. The 
influence of the polytropic index on the shape of 
the distribution can also be seen in the results of 
three-dimensional simulations performed by Li, et 
al. (2003) and Audit & Hennebelle (2010). Stud- 
ies of turbulence in an isothermal medium by Fed- 
errath et al. (2008) indicated that the shape of 
the distribution may also be governed by the char- 
acter of the driving force: solenoidal (divergence- 
free) driving forces result in a symmetrical dis- 
tribution (close to lognormal), while compressive 
(curl-free) driving can create a negative skewness. 
In short, the appearance of our pressure distribu- 
tion seems to favor either 7 > 1 (i.e, somewhere 
between isothermal and adiabatic behavior), tur- 
bulence that is compressive in nature, or some 
combination of the two. 

One important application of our determina- 
tion of the dispersion of thermal pressures is to 
estimate the strength of the turbulence using a 
quantitative comparison based on computer MHD 
simulations. Padoan et al. (1997a, b) introduced a 
scaling relation between the rms dispersion o~ s for 
a log- normal distribution of the quantity s = hip 
in terms of a simple function of the Mach number 
M, 

o- s = [ln(l + 6 2 M 2 )] a5 . (6) 

Investigators that have adopted this formalism 
generally find that their simulations seem to sup- 
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port the validity of a scaling with M shown in 
Eq. [6l but values for the constant b appear to vary 
from one study to the next. Federrath et al. (2008, 
2010) and Brunt (2010) have summarized the out- 
comes for b for many different cases reported in 
the literature: extremes in b have ranged from 0.3 
to 1.0, depending on the conditions in the compu- 
tations. Simulations carried out by Federrath et 
al. (2008, 2010) indicate that whether or not the 
forcing of the turbulence is solenoidal or compres- 
sive can have a strong influence on b. Lemaster & 
Stone (2008) have shown that magnetic fields have 
only a small effect on the relationship between <j s 
and M, 

We can derive a characteristic turbulent Mach 
number for the C I-bearing gas by taking our dis- 
persion for hip, adopting a value for 6, and solv- 
ing for M in Eq. |6l Here, it is appropriate to 
use a volume weighted distribution of pressures, 
since that is the conventional way of describing the 
outcomes of the simulations. A best-fit of a log- 
normal distribution to the portion log(p/fc) > 3 
of the isothermal curve for low I/Io shown in 
Fig. [10] yields a s = 0.46. From the analysis of 
the possible effects of averaging in velocity bins 
that we presented in £)9.31 we acknowledge that 
the true dispersion of log(p/fc) could be as large 
as 0.5, leading to a s = 0.5 In 10 = 1.2. For our 
choice of b, we adopt the finding by Brunt (2010) 
that b = 0.48^0 ii! which was based on the ob- 
served density and velocity variances in cold gas 
with large turbulent Mach numbers in the Tau- 
rus molecular cloud. This value is near the mid- 
dle of the range of those derived from computer 
simulations of MHD turbulence mentioned in the 
above paragraph. With this value for b, we solve 
for M in Eq. [6] and derive M — I.OIqJ for °ur 
lower value of a s and M = 3.7+J'g for the larger 
one. For our representative values /(Ha) = 0.60 
( g3) and T = 80 K (S3), the isothermal sound 
speed c s = 0.50 km s _1 . For the smallest value of 
M minus its error, we expect the velocity disper- 
sion a v = 0.8c s = 0.40 km s" 1 , and for the largest 
M plus its error we expect that a v — 4.8c s = 
2.4km s" 1 . 

Over a wide dynamic range in linear separa- 
tions, the velocity differences for packets of mate- 
rial in the ISM have been observed to scale in pro- 
portion to these separations to a fixed power (Lar- 
son 1979, 1981; Heithausen 1996; Brunt & Heyer 



2002a; Brunt & Kerton 2002). We can factor in 
our values of a v into this relationship to estimate 
the largest characteristic scales for the turbulent 
motions, which in turn indicate the largest cloud 
sizes (or energy injection scales). For the power- 
law relationship, we adopt the recent finding of 
Heyer & Brunt (2004), 

a v = (0.96±0.17)r° c 59±0 07 , (7) 

where r pc is the linear separation in pc. Solving for 
r pc using our velocity dispersions in this equation 
yields r pc = 0.23±g;^| for a v = 0.40 km s" 1 and 
r pc = 4.7^6 for a v = 2.4km s _1 . 

In Fig. [5] we showed the distribution of thick- 
nesses of our C I-bearing clouds for the different 
lines of sight in our survey. The median of all the 
values for the entire collection is 5.5 pc. On the 
one hand, this median value is close to the upper 
end of our range of r pc , which may indicate that 
our larger value of a s , i.e., the largest possible dis- 
persion found in §9.3[ represents the correct value 
for the deviations of thermal pressures. On the 
other hand, the smaller dimensions that apply to 
the direct measurement a s = 0.46 may simply in- 
dicate that we are usually viewing a superposition 
of many independent, smaller clouds along each 
line of sight. We caution that the trend expressed 
in Eq.[7]is evaluated from 12 CO J = 1— emission- 
line data for molecular clouds, which may differ 
from the relationship for the more diffuse regions 
that we have sampled. 

10.2. Time Constants 

Since fluid elements in a turbulent medium have 
physical properties that change with time, it is 
important to establish the time intervals that are 
required for various measurable quantities to con- 
verge nearly to their equilibrium values. There 
are three different contexts where we compare two 
(or more) states of any particular constituent: (1) 
the ratio of C I fine-structure populations, (2) the 
balance between neutral and ionized forms of the 
carbon atoms and (3) the J = to 1 rotational 
temperature Tqi of H2. A fourth time- variable 
quantity of interest is the kinetic temperature of 
the gas, which not only influences the other three 
quantities that we measure but also the manner 
in which the gas responds to disturbances. We 
will compute the characteristic e-folding times for 
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these processes in the following subsections, and 
later we will compare them to the eddy turnover 
times for different size scales. 

In a general context, we can imagine atoms or 
molecules in two possible states with equilibrium 
number densities no, e q in some lower level and 
7ii. oq in an upper one. In equilibrium, 

-Roi^Oxq = -RlO^l.cq ; (8) 

where Rqi and i? 10 are the upward and downward 
conversion rates, respectively. We can propose a 
solution for the time behavior of the lower level, 
no, to take the form 

n (t) = n 0iC q + (no,i - no : cq)e~ 7t ( 9 ) 

as the concentration of no adjusts itself from some 
initial density no,i to its equilibrium value no, eq . 
This form must agree with the condition 

= m(t)R w - n (t)Roi 

= "tot -Rio - n (t)(R w + Rqi) ,(10) 

where n to t — "o (t) + Tlx (t) is the (constant) sum 
of the number densities of the two levels. If we 
insert the proposed time behavior (Eq. [5]) into the 
no(t) term of Eg. 1101 and compare it to an explicit 
differentiation of Eq. |9] with time, we can equate 
the e -7 * terms to reveal that 

"f = R w + R 01 . (11) 

(The sum of the remaining terms without e -7 * 
equals zero.) In essence, any departure from the 
equilibrium level distribution, either positive or 
negative in sign, will decay to the equilibrium con- 
dition in an exponential fashion with an e-folding 
time constant given by Eq. 1111 In the following 
three subsections, we apply this rule to population 
ratios discussed in items (1) to (3) at the begin- 
ning of this section. 

10.2.1. CI Fine Structure Levels 

As a simplification, we consider only the first 
two levels and ignore the existence of the third 
(highest) one. Here, Rqi equals the sum of the 
upward rate constants for various collision part- 
ners times their respective densities. R\q equals 
the sum of the downward rate constants times the 



densities plus also the spontaneous decay proba- 
bility A w . If /l is small, A w = 7.93 x 10" 8 s~ 1 
(Galavis et al. 1997) dominates over the collisional 
excitation (and de-excitation) terms. The inverse 
of A w equals 146 days. If /l is not small, the colli- 
sional terms make Rio + Rqi even larger and thus 
reduce the time constant to less than 146 days. 
Clearly, even for the more complex situation for 
the interactions with the highest of the three fine- 
structure levels, the time constants are extraordi- 
narily short (A21 = 44 days). 

10.2.2. The Photoionization and Recombination 
of Carbon Atoms 

Since the equilibrium concentrations of neutral 
carbon are much smaller than the ionized forms in 
all cases that we consider, as is evident from Fig.Q] 
it is clear that the ionization rate i?oi = {I /IqWq 
dominates over the various recombination terms 
shown in Eq. [T]that make up Riq. Figure [8] shows 
us that I/Iq = 1 is about the lowest value of the 
radiation density that we encounter. Hence, the 
longest time constant that we expect to apply is 
simply Tq 1 = 5 x 10 9 s = 160 yr, and this time 
shortens in proportion to the increase in / above 
the reference value Jo- 

10.2.3. The J = to 1 Rotation Temperature 
T i of H 2 

Cecchi-Pestellini et al. (2005) have performed 
detailed calculations of the time-dependent H2 
level populations in turbulent media that have 
short-lived pockets of hot gas that can leave an 
imprint on the rotation temperatures. Here, we 
focus on a much simpler discussion for Tqi of the 
two lowest rotational levels, since they are relevant 
to our determinations of kinetic temperatures. We 
evaluate the characteristic time for changes in the 
population ratio of J = to that of J = 1 when 
there is a sudden change in the kinetic tempera- 
ture, but we neglect any of the effects of repopu- 
lating the lower levels from cascades from higher 
levels of excitation. 

The rate coefficient for ortho-para transitions 
caused by neutral hydrogen impacts onto H2 is ex- 
tremely low at the temperatures that we consider 
[for T < 300 K, k Q1 < lO^Ws" 1 (Sun & Dal- 
garno 1994)]. For protons, however, the rate con- 
stants are much larger: fcio = 2.0 x 10 _10 cm 3 s _1 



36 




12 3 

Log n(H) 

Fig. 12. — Time constants (in seconds) for (1) the relaxation of the Tbi rotation temperature of H2 to the 
local kinetic temperature (solid contours) , as a result of ortho-para conversions of the lowest two levels due to 
collisions with protons and (2) cooling times t coo i as given in Eq. [12] (dashed lines). The thermal equilibrium 
curve for the CNM that appears in Fig. [T] is shown by the thick curve. This diagram was constructed 
assuming that the gas has /(H 2 ) = 0.6 and He/H=0.09 (see N4.3p . 
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(Gerlich 1990), and k 01 = 9exp(-171/T)fci . The 
solid contours in Fig. [12] show how the time con- 
stants vary over the \og(p/k) vs. logn(H) diagram 
when we combine the rate constants with determi- 
nations of n{p) using Eq. [2] 

10.2.4- The Kinetic Temperature 

Unlike the cases that we considered in § £110.2.11 
110.2.31 for the kinetic temperature we must work 
with a continuous variable instead of a popula- 
tion ratio of two states of some constituent. Thus, 
a somewhat different tactic is needed to assess 
the characteristic relaxation time. Wolfire et al. 

(2003) have evaluated the isobaric cooling time for 
the ISM and find that 

< cool = 7.40xl0 u (—)' ( ' s. 

\80KJ V3000cm- 3 K7 

(12) 

They state that this formula is valid to within a 
factor 1.35 over the temperature range 55 < T < 
8500 K. By itself, the coefficient in Eg. Il"2"l applies 
to conditions very close to our median tempera- 
ture and pressure in the survey (Tqi = 80 K and 
p/k = 3000cm~ 3 K), and it is not much different 
from the relaxation time for Tqi (3.95 x 10 11 s = 
12,500yr) at the same temperature and pressure. 
We depict values of t coo \ by the nearly straight, 
dashed contours in Fig. Q21 

10.3. Turbulent Eddy Crossing Times 

As the length scales become smaller, the dwell 
times for conditions become shorter. It then fol- 
lows that these shortened durations could curtail 
physical stasis in certain respects. To estimate 
the time scales At = r/Av for changes to occur 
in a turbulent eddy with a characteristic radius r, 
we can once again make use of the power-law re- 
lationship between velocity shears Av and length 
scales r pc (as we did in ijl0.1.3[) . but this time 
we adopt the findings taken from observations at 
smaller scales. A slight reduction of the slope 
seems to occur for these shorter length scales: Fal- 
garone et al. (1992) conclude that Av « r pc 4 for 
10~ 2 < r pc < 1, and their result agrees with that 
of Brunt & Heyer (2002b) and Heyer & Brunt 

(2004) at a common scale r pc = 1. This veloc- 
ity trend for the shorter lengths is consistent with 
the theoretical study of turbulence by Boldyrev et 
al. (2002), and we will adopt it for our investiga- 



tion of time scales. We recognize, however, that 
there are isolated observations, such as those car- 
ried out by Sakamoto (2002), Sakamoto & Sunaka 
(2003) and Heithausen (2004, 2006), that show 
some specific regions where the velocity differences 
measured over r pc ~ 10~ 3 — 10 _1 are almost one 
order of magnitude above this velocity-size rela- 
tionship; see Fig. 10 of Falgarone et al. (2009). 
Also, observations of CO emission by Hily-Blant 
et al. (2008) demonstrate that isolated concen- 
trations of turbulent energy over small scales cre- 
ate occasional, large velocity deviations that go 
well beyond the tails of a Gaussian distribution. 
Finally, Shetty et al. (2010) indicate that projec- 
tion effects in the position-position- velocity (PPV) 
data overestimate the power-law slope (by one to 
a few tenths) and underestimate the velocity am- 
plitudes (by about a factor of two) in a 3D phys- 
ical position-position-position (PPP) representa- 
tion of a turbulent medium. With these points 
in mind, we make an extrapolation of the trend 
Av = r pc 4 km s _1 toward very small scales to yield 
At = r pc 6 Myr, but acknowledge that in some 
circumstances this form for At may significantly 
overestimate the time span for rapid changes in 
conditions. 

For a length scale r pc = 0.00046 (or 95 AU), we 
expect that At = 0.01 Myr (i.e., 10 115 s). Along 
the CNM equilibrium curve shown in Fig. [T2j this 
time equals t coo \ at \og(p/k) = 3.85. The crossing 
time is about equal to the e-folding time for Tqi 
at a slightly lower pressure, \og(p/k) = 3.6. Thus, 
in short, for scale sizes smaller than about 100 AU 
(but perhaps a few thousand AU for some of the 
more active regions) we can expect that turbulent 
fluctuations at the pressures that we are consider- 
ing will depart from the CNM thermal equilibrium 
curve (7 « 0.7) and exhibit an effective 7 that is 
somewhere between 0.7 and the adiabatic value of 
5/3 (for a pure atomic gas). Over smaller scales 
(or lower pressures) Tqi may depart from the lo- 
cal kinetic temperature. Over all of the practi- 
cal size scales, the equilibrium results for the C I 
fine-structure excitations and C ionization should 
apply, but with the provision that their outcomes 
depend on the instantaneous temperature. 
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10.4. Possible Effects from Coherent, Large 
Scale Disturbances 

Up to now, we have considered the effects of 
compressions and rarefactions caused by random 
turbulent motions. However, the injection of me- 
chanical energy over macroscopic scales in the 
Galactic disk can also create deviations in pres- 
sure. Supernova blast waves constitute a principal 
source of this energy in the ISM. We know that 
there are strongly elevated pressures inside clouds 
that have recently been overtaken by a supernova 
blast wave, as shown by the enhanced C I fine- 
structure excitations that appear in the spectra of 
stars within and behind the Vela supernova rem- 
nant (Jenkins et al. 1981, 1984, 1998; Jenkins & 
Wallerstein 1995; Wallerstein et al. 1995; Nichols 
& Slavin 2004). We have good reason to expect 
that pressure increases with somewhat smaller am- 
plitudes should persist even within remnants that 
are no longer identifiable because they are so old 
or disrupted. 

In a more general context, at random locations 
in the disk of the Galaxy the outcomes for the 
thermal pressure enhancements arising from the 
effects of supernovae are expected to be appre- 
ciably different for the various broad temperature 
regimes in the ISM, as shown by several differ- 
ent computer simulations (de Avillez & Breitschw- 
erdt 2005a; Mac Low et al. 2005). The simula- 
tions have many free parameters that influence the 
properties of the average pressures and the shapes 
of the distribution functions. For this reason, we 
will restrict our attention to a very simple test 
of the plausibility that, beyond the limited range 
of fluctuations caused by turbulence, there is a 
broader, underlying spread of pressures caused co- 
herent, large scale mechanical disturbances arising 
from supernova explosions. 

We can adopt a tactic similar to one developed 
by Jenkins et al. (1983) in their comparison of 
C I pressures to a prediction based on the the- 
ory of the three-phase ISM advanced by McKee 
& Ostriker (1977). Small, neutral clouds that are 
overtaken by an expanding supernova blast wave 
should rapidly (and adiabatically) adjust their in- 
ternal thermal pressures to equal that of the hot 
medium well inside the remnant's boundary. We 
can now make a simple prediction of what would 
happen if these clouds actually defined the trend 



of the pressure distribution well above the mean 
pressure and then compare this outcome with our 
observations. 

If the radius r of any remnant in the adiabatic 
phase grows in proportion to t v and its volume- 
weighted average internal pressure p is propor- 
tional to r a , we find that 

dp/dt = (dp/dr)(dr/dt) cx r a ~ 1+{r, ' is >^. (13) 

A time-averaged occupation volume V{p) is then 
given by 

V{p) cx r 3 /(dp/dt) = r 3 - Q+1 /", (14) 

which gives an overall volume filling factor per unit 
logp that is proportional to V(p)p = p~ 14 / 9 for 
a = -3 and 77 = 3/5 (McKee & Ostriker 1977), as 
long as the remnants do not overlap each other, 
which should be true at pressures well above the 
median pressure. 

The histogram-style trace in Figure [10] shows 
our thermal pressure distribution on the assump- 
tion that the overtaken clouds contract adiabat- 
ically, i.e., with 7 = 5/3. Unlike the smooth 
curves shown in this figure, this distribution repre- 
sents our entire datasct, i.e., not just the instances 
where I/Iq < 10 5 . Our reason for this choice 
is that we wish to avoid a bias against regions 
of higher than normal starlight intensity, because 
the locations of supernova remnants are correlated 
with those of associations of early-type stars. As 
Fig. llOl shows. the pressure distribution has a slope 
that is roughly consistent with dV/dlnp cx p -14 / 9 . 

10.5. High Pressure Component 

In Sj5]we proposed that some small fraction of all 
of the gas that we observed has an extraordinarily 
high pressure (p/k > 3x 10 5 cm" 3 K, T > 80K),in 
order to nudge the /2 outcomes to locations above 
the normal equilibrium tracks shown in Fig. [2] We 
now explore various explanations for the excesses 
in /2, starting with ones that do not imply the 
presence of small amounts of gas at high pressures. 
Later, on the premise that the existence of the 
high pressure material is indeed real, we review 
some suggestions made by other investigators on 
its possible origin. 

10.5.1. A Misleading Conclusion? 

Before we fully accept our interpretation that 
the anomalously high /2 measurements imply the 
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existence of small amounts of high pressure mate- 
rial, well separated from the main pressure dis- 
tribution function presented in £19. 11 we should 
briefly investigate possible errors in the interpre- 
tation of the outcomes in fl and /2. One possi- 
bility is that the excitation cross sections or the 
decay rates for the excited levels have systematic 
errors that underestimate the populations in the 
3 P 2 state (C I**) relative to those in the 3 Pi level 
(C I*). We feel that this is unlikely, since earlier 
calculations of these quantities that appeared in 
the literature (i.e., the ones adopted by JTOlF^h 
did not yield outcomes that predicted greater val- 
ues of /2 for their respective fl counterparts. 
However, on more fundamental grounds we do not 
feel qualified to comment on the accuracy of the 
atomic physics calculations, so we will not pursue 
this issue further. 

Another possibility for misleading results could 
be errors in our determinations of fl and /2. Two 
possibilities for the origin of such errors could ei- 
ther be errors in the adopted /-values for the C I 
transitions or our under- appreciation of the effects 
of misleading apparent optical depths caused by 
unresolved, saturated substructures in the absorp- 
tion line profiles. For the former of the two, we feel 
that our investigation discussed in Appendix [B] 
provides some assurance that we are not experi- 
encing systematic errors in the relative strengths 
of weak multiplets versus the strong ones. How- 
ever, our derived /-values rely on the correctness 
of the published relative line strengths within mul- 
tiplets. These relative strengths have a direct 
influence on our derived values of N(C I*) and 
N(C I**), relative to each other and to N(C 1). 

As for the possibility that we are being mis- 
led by incorrect optical depth measurements, we 
feel that the precautions that we discuss in ! )A.2.4I 
for screening out such cases provide adequate safe- 
guards. Moreover, it is reassuring to see that 
for individual determinations of the apparent frac- 
tion of C I in the high pressure component, <?high, 
in each velocity bin (i.e., not the overall aver- 
ages shown in Table [3]), there is no trend with 
N(C Itotai), an effect that we would have expected 
to see if the phenomenon were driven by the 
strengths of the absorption lines. 



Appcndix lA.3l discusses our current updates for the atomic 
parameters. 



Still another effect to examine is the possibility 
that there is an anomalous means for exciting the 
fine-structure levels. Positively charged collision 
partners will give proportionally stronger excita- 
tions of the second excited level of C I, as indi- 
cated by the differences in cross sections for pro- 
tons compared to neutrals - see Fig. 1 of Silva 
& Viegas (2002). If ambipolar diffusion (ion- 
neutral slip) created by MHD shocks and Alfven 
waves create enough suprathermal protons (and 
heavy element ions) to further excite the C I, 
they might help to explain the larger outcomes 
for /2. While this is a qualitatively attractive 
explanation, in a quantitative sense it seems to 
fail: the required fractional concentration of the 
positively charged collision partners, greater than 
about 30%, seems to be unreasonably large (e.g., 
the conditions n(H I) = 2 cm" 3 , T(H I) = 600 K, 
n(p) = 0.6 cm" 3 , T(p) = 20, 000 K should give 
fl = 0.23 and /2 = 0.066, which is not far from 
our measured average shown by the white x in 
Fig. [2j Smaller ion fractions fail to do so how- 
ever). 

10.5.2. The Amount of the High Pressure Com- 
ponent 

We now move on to the premise that we ad- 
vocated earlier in Sj5] that the anomalously large 
values of /2 arise from a small admixture of high 
pressure gas in virtually all of the cases that we ex- 
amined. In terms of N(C I), the overall fraction 
5high is usually about 5%. However it is important 
to note that, except in the presence of exception- 
ally strong ionization field strengths, this outcome 
must arise from much smaller proportions of H I 
because the neutral fraction of carbon increases 
with pressure, making small amounts of high pres- 
sure gas far more conspicuous. For instance, we 
can expect a factor 100 enhancement in fractional 
amount of C I, n(C I to tai)/KC II) + n(C I to tai)], 
when the gas is at log(p/fc) = 6, T = 300 K over 
that which would apply to material with more 
conventional physical conditions log(p/fc) = 3.6, 
T = 80 K. On average, this makes the fractional 
mass contribution of H I in the high pressure com- 
ponent only about ^high/100, a 0.05% mass frac- 
tion. We add a caution, however, that this fraction 
could be larger if the actual pressure of the high- 
pressure component is lower than the value stated 
above. 
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10. 5. 3. Radiation from the Excited Levels of Car- 
bon Atoms 

Radiative decay of the upper fine-structure level 
of C II is an important cooling mechanism for the 
ISM. The rate of this cooling per unit mass can be 
monitored by either directly observing the emis- 
sion at 1900 GHz (157.7 ^m) (Stutzki et al. 1988; 
Langer et al. 2010; Pineda et al. 2010; Velusamy 
et al. 2010) or by viewing the C II* absorption 
features at 1037.018 and 1335.708 A, as has been 
done for both the ISM of our Galaxy (Lehner et 
al. 2004) and the distant, damped La systems in 
quasar spectra (Wolfe et al. 2003a, b). It is worth- 
while to ask the question: could the regions that 
we view with enhanced pressures make an impor- 
tant contribution to either the absorption or emis- 
sion measurements? It is difficult to formulate a 
precise answer, since we do not fully understand 
the nature of these regions. On the one hand, the 
factor 100 diminution in the H I content mentioned 
in the previous section is approximately offset by a 
factor 100 enhancement in the collision rate for ex- 
citing the upper level of C II. In this circumstance, 
as long as we are still below the critical density 
for establishing the C II* population, our typical 
value of <7high of 5% would be approximately the 
correct answer for the enhancement of the emis- 
sion intensity (in the optically thin limit) or for 
the increase in N(C II*) over that from the gas at 
normal pressures. On the other hand, if the high 
pressure regions are located at sites where the pho- 
toionizing radiation level is much higher than else- 
where, then the H I concentration is not strongly 
diminished but the population of the upper C II 
levels is still very high. Here, the column densities 
of C II* could be considerably higher than 5% of 
the total and the regions that hold this material 
could emit a substantial amount of radiation. 

It is much easier for us to make a quantita- 
tive assessment of the enhancement of radiation 
from the excited levels of C I because the popu- 
lations of the two upper levels are exactly what 
we observe. The value of fl within the high pres- 
sure gas should be about twice that of the normal 
gas; hence we can expect that the radiation at 
492 GHz (609 fim) seen toward most of the translu- 
cent clouds (Heithausen et al. 2001; Bensch et al. 
2003) should only be increased by about 10%. We 
estimate that the value of /2 in the high pres- 
sure gas is about 15 times as large as that in the 



low pressure gas, so about 44% of the radiation at 
809 GHz (371 /im) could arise from the high pres- 
sure regions. 

10.5.4- Possible Origins of the High Pressure 
Gas 

As we pointed out in <|5j in order to obtain the 
composite fl and /2 values that we found for the 
entire survey, the high pressure component had 
to be a distinct population whose distribution in 
pressure was well removed from the low pressure 
material. We demonstrated in Fig. [3] that it could 
not be a diminishing tail resembling a power law 
that extends away from the main, low pressure dis- 
tribution. In the context of turbulence theory, this 
poses a challenge in the interpretation of the high 
pressure gas, unless one could propose an explana- 
tion for the absence of intermediate mass fractions 
at pressures between the two extremes. 

In order to justify the presence of certain 
molecules in the ISM that require endothermic re- 
actions for their production, such as CH+ , Joulain 
et al. (1998) proposed the existence of hot gas 
concentrations within highly confined dissipation 
regions created by turbulence. From a computer 
simulation, Pety & Falgarone (2000) found that 
extraordinary physical conditions could arise in re- 
gions that were selected to have special dynamical 
conditions, such as larger than normal amounts 
of vorticity or negative divergence. Further stud- 
ies from theoretical or observational perspectives 
have been presented by Godard et al. (2009) and 
Hily-Blant et al. (2008). These investigators con- 
cluded that the volume filling factors for these 
regions are small (a few percent), but not as small 
as the mass fractions that we reported in 8 110. 5.21 
While the intermittent emergence of extreme con- 
ditions within highly confined dissipation regions 
in a turbulent regime is an attractive prospect 
for explaining our high pressure gas, it must nev- 
ertheless satisfy our requirement for a distinct 
separation from the pressure fluctuations arising 
from regular turbulent disturbances instead of a 
continuous, low level extension thereof. 

In fJ5]we showed evidence that the greatest ex- 
tremes in pressure occurred for gases at unusually 
large negative velocities, and this interpretation 
fits in well with the concept that the target stars 
(and their neighboring stars) play a role. Indeed, 
an inspection of Table [3] shows that in some di- 
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rections, several adjacent sight lines all show ele- 
vated pressures compared to the typical pressures 
derived from the full sample. Two prominent ex- 
amples are stars near or within the Carina and 
Orion Nebulae. These are dynamically disturbed 
regions, and they are also regions of significantly 
enhanced starlight density. This supports the no- 
tion that the stars somehow raise the pressures in 
their surroundings and that high values of I /Iq in- 
dicate both recent, enhanced star formation and a 
more highly pressurized ISM. 

Figure [T3l indicates that the quantities ghigh and 
I/Io also seem to be connected to each other. Gen- 
erally, we can see that cases where <?high > 0.2 
appear to require that log(///o) > 0.5 and that 
there were very few outcomes that had <?high < 
0.2 that had log(//io) > 1-5- It is unclear 
whether the dominant cause for pressurization 
is from interactions with mass-loss ejecta, small 
clouds experiencing a photoevaporation "rocket 
effect" (Bertoldi 1989; Bertoldi & McKee 1990; 
Bertoldi & Jenkins 1992), or the sudden creation 
of an H II region, all of which can compress the 
ambient material and accelerate it toward us. An 
additional possibility is that H I gas near the stars 
is heated more strongly by the photoelectric ef- 
fect from grains (Weingartner & Draine 2001b), 
which could cause a short-term spike in pressure. 
While these effects (or combinations thereof) may 
be the dominant cause for the isolated cases that 
show large values for (/high, we still find significant 
amounts of high pressure material at large posi- 
tive velocities and even small admixtures of high 
pressure material at all velocities. These outcomes 
indicate that other mechanisms unrelated to the 
target stars may play role as well. 

Field et al. (2009) have proposed that the re- 
coil of H atoms following the photodissociation of 
H 2 at the edge of a molecular cloud can create an 
external pressure that helps to confine the cloud. 
They estimated that at locations where ambient 
UV field intensity reaches I/Iq « 60 the recoil 
pressure can be approximately 1.3 x 10 5 cm~ 3 K. 
(Our I / Iq is approximately the same as their % 
intensity parameter.) If this mechanism could in- 
crease the thermal pressure in the outer portions 
of translucent clouds that still have reasonable 
concentrations of H2 and explain the existence of 
the high pressure component that we observe, we 
would indeed expect to see a positive relationship 



between <?high and I /Iq. 

Even though the points in Fig.[13]at first glance 
seem to favor the recoil hypothesis as a possible 
explanation for the high pressure component, our 
enthusiasm for accepting this proposition must be 
tempered by two considerations: (1) our fiducial 
pressure log(p/fc) > 5.5 (a lower limit which is to 
some degree arbitrary and might be relaxed to a 
slightly lower level) requires a value for I /Iq con- 
siderably greater than 60, and (2) the correlation 
seen in Fig. [13] may be a byproduct of some other 
physical effect that responds to higher than nor- 
mal intensities and generates the high pressures. 
For instance, we know from Fig. [8] that values of 
log(p/fc) in the low pressure regime are likewise 
correlated with the intensity of starlight, and the 
greater influence of intermittent dissipation effects 
in the more strongly driven turbulence may ac- 
count for the more conspicuous presence of high 
pressure gas. 

11. Summary 

We have presented a comprehensive analysis 
of HST UV spectra stored in the MAST archive 
for 89 stars that were observed with the E140H 
mode of STIS, with the goal of measuring the pop- 
ulations of the three fine-structure levels of the 
ground electronic state of neutral carbon atoms in 
the ISM. The ultimate purpose of these measure- 
ments was to synthesize a distribution function for 
the thermal pressures in gases that mostly repre- 
sent the cold neutral medium (CNM) in our part 
of the Galactic disk. This work builds upon a sim- 
ilar study of 21 stars in restricted portions of the 
sky carried out by JT01 in a special observing pro- 
gram dedicated to this purpose. We have repeated 
the basic analysis protocol developed by JT01 for 
unraveling the velocity profiles for carbon atoms in 
the separate fine-structure levels from the blended 
features in many different multiplcts, but with a 
few improvements in methodology that are out- 
lined in various subsections within Appendix |A1 

We interpret most of the variations in the out- 
comes for thermal pressures to arise from fluctua- 
tions caused by interstellar turbulence. Some ad- 
ditional pressure excursions that are large in mag- 
nitude but for small mass fractions probably arise 
from the the random passages of infrequent but 
strong shocks created by either stellar mass loss 
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Fig. 13. — The relationship between the fractional amount of high pressure gas, (/high, and the starlight 
intensity relative to the Galactic average, I flo- 



or supernova explosions. 

The following conclusions have emerged from 
our study of C I fine-structure excitations: 

1. The relative populations of the two excited 
fine-structure levels are influenced in differ- 
ent ways by the local physical conditions, 
since the levels have significantly different 
collisional rate constants and energies. This 
feature allows us to sense in any one veloc- 
ity channel the presence of admixtures of 
gas that have markedly different conditions. 
While there is a multitude of possibilities 
for explaining any particular combination of 
level populations, we find that when the data 
are viewed as a whole, the most straight- 
forward interpretation is that practically all 
of the gas in the normal range of pressures 
(10 3 < p/k < 10 4 cm~ 3 K) is accompanied 
by very small amounts (of order 0.05%) of 
gas at anomalously large pressures and tem- 
peratures (p/k > 10 5 - 5 cm~ 3 , T > 80 K). In 
a small fraction of cases, the proportion of 
the gas at high pressures is markedly higher 



than this level, both because the amount of 
C I is greater and the local radiation den- 
sity is high (which makes more of the carbon 
atoms singly ionized). 

2. For a substantial number of our lines of 
sight, we can make use of molecular hy- 
drogen rotation temperatures Tqi between 
J = and 1 to define the local kinetic 
temperature. Such temperatures are use- 
ful in defining one of the free parameters in 
solutions for the level populations. As an 
added benefit, we can explore whether or 
not, in a general statistical sense, the pres- 
sure outcomes are related in some way with 
such temperatures. We find only a weak an- 
ticorrclation, which indicates that pressure 
fluctuations do not appear to be the domi- 
nant cause for temperature changes from one 
place to the next. 

3. Excluding the small amounts of high-pressure 
gas mentioned earlier, the pressures of most 
of the CNM material show some correlation 
with the local radiation densities, as sensed 
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by the observed ratio of O I (or sometimes 
S II) to C I followed by an application of 
the equation of ionization equilibrium with 
plausible values for O/C and S/C to derive 
N(C II). We interpret this trend as arising 
from the fact that the stars that create this 
radiation are sources of enough mechanical 
energy to make the pressures higher than 
normal. 

4. The main part of the mass- weighted distri- 
bution of pressures in our complete sample 
approximately follows a lognormal distribu- 
tion with a mean value for log(p/fc) equal to 
3.58 and a standard deviation of 0.175 dex. 
However, for log(p/fc) < 3.2 or > 4.0 the 
amount of material is greater than a contin- 
uation of the lognormal distribution. 

5. In order to sense the distribution of pressures 
in regions well removed from the sources of 
mechanical disturbance (i.e., the stars that 
emit large amounts of radiation), we have 
isolated for study only those cases where the 
radiation density is less than 10 05 times the 
overall average level. Under these circum- 
stances, the tail on the high pressure side 
of the distribution becomes suppressed, and 
the remaining distribution develops a nega- 
tive skewness. We supply a polynomial ex- 
pression that fits this distribution (Eq. [4]) 
that is shown in Fig. [9] About 23% of the 
material in this distribution is below the 
minimum pressure for the thermal equilib- 
rium curve of a static CNM in our part of the 
Galaxy, suggesting that short-term fluctua- 
tions in pressure can occur without the gas 
being transformed to a stable warm neutral 
medium (WNM). 

6. The thicknesses of the regions that we were 
able to probe, as measured by the hydro- 
gen column densities divided by their space 
densities, are generally less than 20 pc. The 
filling fractions for the sightlines are gener- 
ally less than 1%. The remaining 99% of a 
typical sightline is filled with much hotter 
gas having densities that are far too low to 
create enough C I for us to measure. 

7. We recognize that even with the over- 
determination of conditions provided by the 



two fine-structure levels, we can still under- 
estimate the dispersion of pressures because 
we are viewing at each velocity an aver- 
age pressure for the superposition of regions 
that could have vastly different pressures. 
We have studied how the dispersions scale 
in proportion to 7V(C It. o tai)~ ' 5 and find 
that the ISM could conceivably be composed 
of independent packets of gas with a true 
rms dispersion in log(p/fc) that could be as 
large as 0.5 dex, which is considerably wider 
than the distribution that we constructed di- 
rectly from the data. The characteristic col- 
umn density of each packet would be about 
N(C Itotai) = 2 x 10 12 cm~ 2 . However, an al- 
ternate interpretation, and one that is quite 
plausible, is that small volumes of gas have 
intrinsic pressure variances that are larger 
than for coherent, larger volumes that might 
be more resistant to perturbations from tur- 
bulent forces. This phenomenon could con- 
ceivably produce the same linear scaling of 
pressure dispersions against N(C Itotai) -0 ' 5 
that we observed. 

8. On the basis of our findings reported in items 
4 and 7 above, we derive characteristic tur- 
bulent Mach numbers for the C I-bearing 
gas to range between 0.8 and 4.8. Since 
the speed of sound is about 0.5 km s -1 if 
T = 80 K, we expect that the 3-dimensional 
velocity dispersion <j v to range between 0.40 
and 2.4 km s _1 . If we equate these numbers 
to observations of velocity structure func- 
tions in the ISM, we find that the charac- 
teristic size r of the clouds or the outer driv- 
ing scale of the turbulence is probably in the 
range of approximately 0.2 < r < 4.7 pc. 

9. Gas with radial velocities well outside the 
range of motions expected for differential 
galactic rotation is more likely than usual to 
exhibit exceptionally large pressures. This 
link of pressures with kinematics helps to 
support the interpretation that shocks and 
turbulence play an important role in creating 
the positive excursions in pressure. Pack- 
ets of C I moving at negative velocities show 
larger pressure excursions than for those at 
positive velocities. We explain this differ- 
ence in terms of an observational bias that 
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favors our viewing the near sides of pressur- 
ized shells that are expanding away from our 
target stars. 

10. There is a broad range of time scales that are 
needed to reach equilibrium values for var- 
ious quantities and physical processes that 
are relevant to our study. From the short- 
est to the longest they are as follows: (1) 
C I fine-structure level populations (of or- 
der 100 days), (2) the balance between C I 
and C II established by the competition be- 
tween photoionizations and various means of 
recombination (160 yr, or shorter if the radi- 
ation density is larger than average), (3) the 
coupling of the J = to 1 rotation temper- 
ature of H 2 to the local kinetic temperature 
(10 4 yr for typical conditions: \og(p/k) = 3.5 
and T — 80 K), and (4) the cooling time 
for the ISM (3 x 10 4 yr for the same condi- 
tions) . We compute the eddy turnover times 
for turbulent eddies having a radius r using 
a relation At = r/Av with an extrapola- 
tion to small scales Av = Tp C 4 km s _1 , and we 
find that the only items of consequence for r 
smaller than about 100—1000 AU are (3) and 
(4). Over these extremely small scales, de- 
lays in the adjustments of H2 rotation tem- 
peratures will give misleading readings for 
the local kinetic temperatures, but the dif- 
ferences in the two temperatures should be 
minor, especially since we can use T i only to 
indicate an average temperature over many 
small volumes. Likewise, any lag in the 
thermal response of the gas will make its 
polytropic index 7 closer to the adiabatic 
value, rather than matching the slope of 
the thermal equilibrium curve for the CNM 
(7 w 0.7). The fact that this may be hap- 
pening is supported by the negative skew- 
ness of our distribution in log(p/fc), which 
indicates that the turbulent fluctuations are 
consistent with 7 > 1. 

11. For logp/k above 4.0, we find a slope in 
the relationship between the logarithms of 
the volume fractions of the gas and \og(p/k) 
to be consistent with a power-law slope of 
— 14/9 that is expected for random penetra- 
tions of expanding supernova remnants in 
various stages of development. 
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A. Improvements in the Analysis over that of JT01 

Over the time since the publication of our initial survey of thermal pressures (JT01), we have had the 
benefit of some extra opportunities to recognize various ways to improve our data reduction and analysis 
methods. In the following subsections, we discuss these new features in our more refined treatments, and we 
also cover a number of improvements in some basic atomic data that have emerged since the earlier study. 
We refrain from discussing here all of the fundamentals of how the analysis was carried out; JT01 explained 
this in some detail. 

A.l. Available Transitions and Minor Modifications of Some /-Values 

The observing program of 21 stars carried out by JT01 had a wavelength coverage that included all of the 
multiplets out of the ground electronic state of C I at wavelengths longward of 1188 A, except for Multiple!^ 
3 centered at 1561 A. Their analysis was based on intensity profiles recorded for all of these multiplets, except 
for a few cases where either the lines were too weak to be useful or suffered interference from other atomic 
species (see Table 3 of JT01). They adopted a revised set of /-values for the transitions that gave internally 
consistent outcomes for their analysis, based on comparisons of optical depths described in §5.3 of their 
paper. 

Some of the observations used in the current study have exceptionally high signal-to-noise ratios. We have 
made use of these results to further revise (or determine for the first time) the /-values for some transitions 
that were too weak to evaluate accurately in the earlier study, again by implementing the analysis of JT01. 
We list in Table [5] our newly adopted /-values for the current study For lines not listed in this table, we 
used the values given by JT01. New /-values were also determined for Multiplet 3, which was covered by 
observations in the archived data but not in our original survey. New /-values determined here are very 
similar to those derived by JT01; both sets differ appreciably from ones published in the literature. In 
Appendix [B] we describe a special study that supports the validity of both our new and older /-values. 

A. 2. Improved Quantitative Estimates of the Uncertainties 

A critical aspect of our study of the thermal pressures is the proper understanding and control of errors, 
since, if they are large enough, they can mislead us into thinking the distribution function is broader than 
reality. For this reason, we have instituted a number of improvements in developing quantitative estimates 
for various sources of error so that we can more reliably screen out measurements of inferior quality and 
have a better confidence that the remaining errors are inconsequential. The discussions in the following 
subsections build upon the concepts presented by JT01 (see their Sections 4 and 5). 

A. 2.1. Relative Errors in the Observed Optical Depths 

For absorption features that are moderately or very strong, the principal source of uncertainty in any 
intensity measurement is that produced by random noise fluctuations in the counts of photoevents. In the 
study of JT01, the effects of these errors were propagated through the analysis, serving as a guide on the 
most appropriate weight factors for intensities recorded in different multiplets. JT01 also recognized the 
existence of uncertainties in the adopted background level by using a formula (their Eq. 7) that reduced or 
eliminated the relative weights of stronger features that came close to zero intensity. For very weak lines, 
however, systematic uncertainties in the establishment of the continuum level are also important. 

In the current study, we now add the continuum uncertainties to the noise errors in quadrature. This has 
practically no effect for moderately strong features, but it now decreases the relative importance of very weak 
multiplets in the final solutions. As did JT01, we adopted best fitting Legendre polynomials for the continua 
at locations somewhat removed from the features. To construct the probable errors in these continua, we 
evaluate the expected errors in the polynomial coefficients, as described by Sembach & Savage (1992), and 



Multiplet numbering system from Moore (1970), also adopted by Morton (2003). 
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Table 5 
New /-values' 



Multiplet b 


A 

(A) 


log (/A) 


2 


1656.267 


2.012 c 




1656.928 


2.392 c 




1657.008 


2.266 c 




1657.379 


1.788 c 




1657.907 


1.914 c 




1658.121 


1.789 c 


3 


1560.309 


2.312 




1560.682 


2.187 




1560.709 


1.710 




1561.340 


1.488 




1561.367 


0.311 




1561.438 


2.236 


7.01 


1277.190 


-0.021 


12 


1192.451 


0.711 




1192.835 


0.129 


14 


1189.447 


1.006 




1189.631 


1.246 



a Changcs from values adopted by 
JT01 

b Multiplet numbering system 
from Moore (1970), also adopted by 
Morton (2003). 

c Adopted from Morton (2003). 
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then we multiply them by two in order to make an approximate allowance for additional uncertainties caused 
by the arbitrariness in selecting the most appropriate polynomial order. From the sizes of residual errors 
that seem to be extended over broad ranges of velocities where no absorption is evident, this global increase 
in the estimates for the continuum errors seems to be appropriate. 

A. 2.2. Error Estimates for N (C I), N (C I*), and N (C I** ) 

In their Section 5, JT01 discussed the various sources of both random and systematic errors in column 
densities. They estimated the combined magnitudes from most of these errors by measuring the amplitudes of 
fluctuations in N at velocities well removed from obvious C I features. These empirical determinations should 
be satisfactory for velocities where the amplitudes of the absorptions are small, but they underestimate the 
errors at locations where the intensity levels are well below the continuum. In these circumstances, noise 
deviations can create much larger uncertainties in the combined optical depths for different lines within a 
multiplet, and these enhanced uncertainties propagate their way through to the solutions for the column 
densities. 

We have now implemented evaluations of errors that should apply equally well to both the weak and strong 
portions of the absorptions. Eq. ( q of JT01 was used to construct a design matrix from which least-squares 
solutions would emerge for the three column densities at each of the different velocities. The inverse of this 
matrix gives the expected variances and covariances of these variables, as long as there is a proper scaling of 
the rms errors in all of the optical depths that contribute to the matrix terms. By incorporating the error 
derivations discussed in the above section f ^A.2.11) . we believe that our evaluations of the uncertainties in 
the optical depths <J T u) properly include all effects except for some that are unquantifiable, such as possible 
detector artifacts, distortions in optical depths caused by unresolved saturated structures (discussed in EIA.2.41 
below), and errors in the adopted /-values of the transitions. The magnitudes of the covariances are much 
smaller than the variances, so the squares of the errors can be extracted simply from the terms in the main 
diagonal of the inverse matrix. 

In order to obtain an approximate validation that the error calculations give reasonable results, we can 

draw upon three observations of the star HD2I9I88 taken at different epochs. If we exclude the velocity 

component that has been identified by Welty (2007) to vary with time, we find that for the expected errors 

that were calculated according to the prescription given above, the dispersion about the weighted mean 

of individual observations of column densities at different velocities yielded a x 2 = 750 for 536 degrees of 

freedom. This evaluation was performed after the outcomes with a spacing of 0.5 km s _1 had been binned 

by a factor of 3, so that the samples were approximately commensurate with the velocity widths of the 

detector's pixels. One can surmise from this study that our error calculations probably underestimate the 

l 

true errors by a factor of approximately (536/750)2 = 0.85. There were no obvious differences between the 
magnitudes of individual \ 2 outcomes within strong C I absorptions as opposed to those outside the features. 

Over velocity intervals where we sense that there is no C I absorption, we detect some very low level, 
smooth deviations that are still present and that add to the random short-scale noise caused by photon 
counting statistical variations. These deviations are probably caused by slight inadequacies in the fitting of 
Legendre polynomials to the true continuum levels, which can exhibit troublesome variations for stars with 
low projected rotational velocities. These deviations have a relatively minor influence on column densities at 
specific velocities, but they can contribute nonnegligible errors in iV(C Itotai) integrated over all velocities, 
since these errors are coherent from one velocity to the next and thus can build up in a linear fashion. In 
computing the errors associated with the values of N(C Itotai) listed in Table [3] we took such errors into 
account (by direct addition across velocity channels rather than quadrature sums). We also allowed for the 
fact that the smaller scale statistical errors were coherent over contiguous stretches of 2.6 velocity channels, 
which represent the width of a pixel on the STIS ultraviolet detector 



The constant A in this equation should actually be C, where C was denned on the preceding page (and used in Eq. 3). 
This channel width should not be confused with the line spread function of the spectrograph, which is about twice as large. 
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A.2.3. Errors in fl and /2 



JT01 used their own judgment in making the choices for velocity ranges over which the C I absorptions 
were deemed to be strong enough to include in the presentations of fl and /2. The existence of a few errant 
fl and /2 points in Fig. 7 of JT01 shows that this selection was not ideal. We have thus introduced a more 
rigorous and quantifiable selection criterion for data which are to be regarded as acceptable. 

If we re-express the column densities for C I in the three fine-structure levels as Nq = N(C I), Ni = 
N(C I*), and N 2 = N(C I**), along with a designation for the total column density TVtot = Nq + Ni + N 2l 
then the uncertainties <r(7Vo), cr(TVi), o~(N 2 ) contribute to an overall error in fl in the following manner: 



Likewise, 



a(N ) 



dfl 



dNo 



a(iVi) 



dfl 



a(N 2 ) 



dfl 

dN 2 




{aiN,) 2 [JVtat - N,} 2 + N 2 [a(N ) 2 + a(N 2 ) 2 ] } 



a(f2) = N t - 2 [<t(N 2 ) 2 [N tot - N 2 } 2 + Ni [o-(N ) 2 + a^) 2 ] } 



(Al) 



(A2) 



For a determination of fl and /2 at any particular velocity, we insisted that in order to be considered, 
both of their ler errors had to be less than 0.03. 



A. 2. 4- Sensing Possible Distortions Caused by Unresolved Saturated Portions of the Absorption Profiles 

Even though the spectra considered here had a minimum wavelength resolving power R — A/ A A = 114, 000 
for the E140H grating with the standard entrance sliiF^I (Proffitt et al. 2010), there remains a possibility that 
for some cases we will encounter collections of unresolved, saturated absorption features. In such situations, 
the apparent optical depths will underestimate the true optical depths after instrumental smoothing. One 
can normally detect this condition by noting that strong lines show smaller column densities than weaker 
ones and then applying a correction scheme proposed by Jenkins (1996). Here, however, the correction is 
no longer straightforward because the individual apparent optical depths lose their identity after they have 
undergone the transformations that are needed to unravel the overlaps of C I, CI*, and C I** absorptions 
in each multiplet. 

While we were unable to perform corrections to restore the apparent optical depths to their true (but 
smoothed) values, we could nevertheless detect circumstances where the representations are likely to be 
inaccurate. The implementation of Eq. 7 of JT01 serves to limit the influence of stronger lines that have 
low residual intensities. By varying the threshold intensity parameter I t in this equation, we can shift the 
response of the solutions for the column densities either toward or away from the stronger lines. In so doing, 
if we find that the column density solutions change, we can surmise the outcomes are not stable and that 
distortions are indeed happening. For adopting or rejecting a set of column densities at any given velocity, 
we adopted the following test: If any column density derived using twice the standard value for I t (see JT01 
§5.2.1) came out to be more than 1.2 times that derived from the standard It, we considered the distortion 
to be unacceptably large and the result for the particular velocity was rejected. Otherwise, deviations less 
than the factor of 1.2 were deemed to be acceptable. Stars for which the optical depth distortions seemed 
to be evident are identified in Table [3] (see note a of the table for details) . 

A. 3. New Atomic Physics Parameters 

With the passage of time, re-evaluations of atomic parameters are carried out (and we presume that they 
are better than the older ones). Since the time of publication of our earlier results (JT01), new reaction rates 



R = 200, 000 for the narrow entrance slit and half-pixel intensity sampling used for the stars observed by JT01. 
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for the excitation of the excited fine-structure levels of C I and O I by atomic hydrogen have been published 
by Abrahamsson et al. (2007). We have incorporated these new rates into our interpretations of densities 
and temperatures from the C I level populations. (The O I excititations are discussed in Appendix 151 ) 
Rates for other collision partners are the same as those adopted by JT01. Also, we have now replaced the 
old spontaneous radiative decay rates for the excited levels of the two neutral species by those given by 
Galavis et al. (1997). 

The effect of these changes is that for the representative conditions \og(p/k) — 3.6 and T — 80 K the 
derived values of the thermal pressure increase by 0.05 dex above what we would have obtained from the 
older numbers used by JT01. 

A. 4. Revised Optical Pumping Rates 

As pointed out by de Boer & Morton (1974), the ultraviolet transitions from the C I electronic ground 
state that can be used to measure the column densities of this atom also act to populate the upper fine- 
structure levels through optical pumping by the ambient starlight radiation field. While this effect is usually 
small compared to the excitations by collisions, it nevertheless is a process that should not be ignored since 
there are occasions when the starlight field can be found to be considerably above average. Jenkins & Shaya 
(1979) calculated the rates of optical pumping for the average level of starlight radiation computed by Jura 
(1975a, b) and Witt & Johnson (1973). These rates were used by JT01. 

Our current analysis of the effects of starlight presented in $4.1\ adopts as a standard the more recent 
average radiation intensity Iq(X) specified by Mathis et al. (1983). Here, we used our ionization equilibrium 
calculations for carbon atoms to determine the local starlight intensities in terms of a multiplier I / Iq times 
this field strength. In order to make our analysis consistent with the new standard intensity, and also to be 
consistent with our new /-values for the C I transitions, we have recomputed the matrix for the pumping 
rates between the different fine-structure levels. This new matrix, presented in Table [5] [cf. Table 3 of 
Jenkins & Shaya (1979)], expresses the transition rates expected for the standard field density. For every 
measurement of C I and its fine-structure populations, we adjusted the pumping rates in proportion to the 
corresponding value of I/Iq - see SJH 

Our application of the new pumping matrix overlooks two effects that can modify the the pumping rates 
in opposite directions. First, we have ignored weak C I transitions at wavelengths shorter than 1188 A. An 
approximate compensation for this is our neglect of the self shielding of the stronger transitions deep inside 
some of the most strongly absorbing clouds. 

The overall effect in implementing the revised pumping rates is to raise the median pressure by about 
0.05 dex for most of the C I data. For measurements that apply to cases where I/Iq < 10 5 , as shown in 
Fig. IH1 the change is almost negligible. 

B. A Validation of our Previous Determinations of the C I /-Values 

One aspect of our earlier study (JT01) of interstellar C I that is highly relevant to the work done here 
was our re-evaluation of the relative /-values of the C I lines (for all three fine-structure levels) from one 
multiplet to the next. Initially, we found that the /-values given in the literature did not yield self consistent 
outcomes for the predicted line strengths, after we had derived preliminary values of the C I, C I* and C I** 
column densities. As a result, we derived a new set of /-values that gave the best internal consistencies for 
the strengths of the lines for all of the targets collectively. These relative /-values were based on the 1656 A 
multiplet as a calibration, in our belief that theoretical calculations of the strength of lines in this strong 
multiplet were probably the most reliable. 
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Table 6 

New Optical Pumping Rates Between the C I Fine-Structure Levels (l(n 10 s~ 1 ) a 



Initial J 




Final J 







1 


2 







5.06 


3.99 


1 


1.69 




5.23 


2 


0.08 


3.14 





a Computcd for an optically 
thin medium in the presence of 
an average ultraviolet radiation 
field density in the ISM at a 
Galactocentric radius of 10 kpc, 
as specified by Matins ct al. 
(1983). 
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Our revised /-values showed a steady divergence from the earlier, published ones as the lines became 
successively weaker. For each —1 dex change in a published line strength log/A pu b., the discrepancy 
log/AjToi — log/Apub. increased by about +0.3 dex (see Fig. 3 of JT01). More recent determinations 
of C I /-values have been compiled by Froese Fischer & Tachiev (2004) (henceforth FT04). A plot of these 
newer values against those derived by JT01 still shows the divergence seen in Fig. 3 of JT01, but with less 
scatter in the individual points representing the weakest lines. Henceforth, we will adopt the values of FT04 
as a proxy for /-values determined by investigators other than JT01. 

In our earlier study, we recognized that the sense of the divergence was consistent with the possibility that 
we could have been misled by the underrepresentations of smoothed optical depths for the most saturated 
portions of the strong lines (while at the same time, optical depths of the weak absorptions would be 
properly measured). Such an effect can occur when the lines have very narrow component structures that 
are saturated and unresolved by the spectrograph (Savage & Sembach 1991; Jenkins 1996), and indeed this 
seems to happen in some cases - see note a in Table [3] 

To overcome our worry about these possible distortions in the apparent optical depths, we carried out 
a number of investigations that indicated that they were unlikely to be the explanation for the divergence 
that we found (see the discussion in Section 5.3.2 of JT01). Nevertheless, it is still troubling that theoretical 
and experimental derivations of the C I /-values by different investigators, such as those cited by Morton 
(Morton 2003), consistently gave results that disagreed with the values derived by JT01. Recent evidence 
that interstellar absorption features arising from neutral atoms can occasionally exhibit extraordinarily low 
turbulent velocity dispersions (Dunkin & Crawford 1999; Price et al. 2001; Knauth et al. 2003; Meyer et al. 
2006) reinforces our concerns about possible misrepresentations of averaged true optical depths. If nearly 
all of the C I absorptions that we observed consisted of clusters of needle-like features that had velocity 
dispersions b < 0.4 km s^ 1 created by low temperatures and low turbulent (or shear) velocities, and they 
are separated 17 ! by more than about 2.56, we could have in principle erroneously introduced systematic 
errors that could explain the observed differences between our /-values and the published ones, even though 
the lines were observed at a wavelength resolving power R — 200,000 (see Appendix A of JT01). This 
phenomenon would have had to operate in a consistent fashion from one case to the next, since we almost 
always found reconstructed profiles (using our /-values) to show good agreement with the observations of 
both strong and weak lines. The uniform persistence of this misrepresentation effect seems rather unlikely, 
but we have no grounds for saying that it is completely impossible. 

We now present a single example that disproves the proposition that our results were incorrect because 
we were continually misled by narrow substructures in the profiles. The substance of our argument is that 
for one special case we can compare the observed and predicted absorption features from a velocity complex 
for which the values of b caused by thermal Doppler broadening for any subcomponents must not be very 
much less than the resolution of the STIS spectrograph. 

A cluster of absorption profiles with a peak at about —34 km s _1 toward HD 210839 (A Cep) is unusual 
because it shows exceptionally strong features from O I* and O I**, as shown in Fig. Qj] and the C I 
fine-structure excitation indicates that it is at an unusually high pressure. The ratio n(0 l*)/n(0 I**) 
is a good indicator of the local kinetic temperature of the gas. Figure [TS] shows how this ratio should 
vary with temperature. This relationship follows from solutions for the equilibrium equation based on the 
H I collision rate constants as a function of temperature given by Abrahamsson et al. (20070 and the 
radiative decay rates of Galavis et al. (1997). We find from the absorptions produced by the 1304.86 A and 
1306.03 A transitions that over the velocity interval -40 < v < -20km s" 1 , log iV(0 I*) = 13.410 ± 0.018 
and log AT (O I**) = 13.359 ± 0.023. The nominal value of 1.12 for the ratio of the two is shown by the 
lower of the two horizontal dashed lines in Fig. [15j while the other dashed line just above it indicates our 



7 Narrow features that have small separations and thus are partly blended are of no concern, since they behave much like features 
with much larger b values. 

8 Abrahamsson et al. did not list rates for T > 1000 K. The collision rate constants that we needed for T somewhat above 1000 K 
are simply extrapolations from their values just below 1000 K. 
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estimate of 1.35 for an upper limit based on the +2cr error for N(0 I*) and — 2a limit for N(0 I**). It 
follows that our favored value for T is 390 K, and our most conservative (i.e., low) value is 200 K, both of 
which arc substantially lower than an upper limit T < 660 K based on a fit that gave b = 0.53 km s" 1 for 
the strongest component of the K I line in this same velocity complex observed at a resolving power of 
0.56 km s _1 (FWHM) by Welty & Hobbs (2001). (The weaker K I component at a more positive velocity 
yielded a b value that indicated that T < 1550 K.) Since the thermal Doppler broadening for T — 390 K is 
expected to yield a value for b of only 0.41 km s _1 , there is either some extra broadening due to turbulence or 
the profiles might be split into even smaller separate components that were not resolved by Welty & Hobbs. 

One might argue that the near agreement of the O I* and O I** column densities could arise from 
these features themselves being strongly influenced by unresolved, saturated absorption spikes within the 
velocity complex. Extreme saturations of this sort would tend to deceive us into deriving nearly the same 
column densities for the two species. However, we can test for this by examining weaker features from 
the 1040 A multiplet covered by the Far Ultraviolet Spectroscopic Explorer (FUSE). One such spectrum of 
HD 210839 is available in the MAST archive. No absorption by the O I* 1040.94 A transition is readily 
apparent in the spectrum; a formal measurement of the equivalent width over the appropriate velocity range 
yields 0.8 ± 1.3 mA. For the column density of O I* that we derived, one would have expected to find 
VF A (1040.94 A) = 2.2± 0.1 mA. If we had underestimated N(0 T) because the 1304.86 A feature had hidden 
saturation, the equivalent width of the weaker feature in the FUSE spectrum would have been larger than 
this expectation. We are unable to perform the same test for O I** because the 1041.69 A feature suffers 
from interference from the Lyman 6—0 R(6) transition of H2 at nearly the same wavelength. 

A separate argument that disfavors a strong internal saturation of the 1304.86 A feature and a less strong 
effect with the 1306.03 A feature is that the shapes of the two absorptions do not differ from each other 
appreciably. If anything, the O I* absorption seems more strongly peaked in the center, an effect opposite 
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Fig. 15. — The expected ratio of n(0 l*)/n(0 I**) for collisions in an environment where n(H I) < 100 cm~ 3 
(solid curve) and n(H I) = 1000 cm -3 (dashed curve). For the velocity component at —34km s _1 in front 
of A Cep, our nominal value and upper limit for N(0 T)/N(0 I**) are indicated by the horizontal dashed 
lines, and the temperatures T — 200 and 390 K that apply to these values are shown by the vertical dashed 
lines. 



to an expectation that saturation could be occurring in the strongest part of the profile. This effect might 
indicate the presence of slightly cooler gas at velocities very near the central velocity of —34 km s _1 (the 
ratio of apparent optical depths at this central point is 1.42, which would indicate T ~ 120 K if we were to 
use the curve shown in Fig. I15|) . 

Returning to the topic of subcomponents within the C I, C I* and C I** profiles, for our conservative 
lower limit T = 200 K we expect the Doppler broadening to produce a b value for carbon atoms of at least 
0.53km s~\ which is somewhat less than the STIS instrumental profile function with b = 0.90km s" 1 (if 
its shape is Gaussian). We can build a model of this complex composed of separate components, each with 
b = 0.53 km s _1 , on top of a low amplitude, broad shoulder. This model has been tailored such that, after 
it has been smoothed by the STIS instrumental spread function, it duplicates exactly the shape of the C I* 
profile that we derived. It also creates structures that are consistent with the 3 K I component parameters 
derived by Welty & Hobbs (2001), after we acknowledge the fact that the thermal broadening of the K I 
should be less than that for C I. 

If we could have resolved the partly blended complex of profiles perfectly, we would have found that the 
apparent optical depth r a (w) [in this case equal to the true optical depth t(v)] divided by /A of the transition 
should remain constant for different transitions at all velocities v, regardless of line strength. However, after 
the profile function has been smoothed by the instrumental profile, there can in principle be some violation 
of this equality. Our model was constructed in a manner to create the worst conceivable situation that would 
aggravate this violation, i.e., we adopted the lower limit for T instead of the nominal one, and we assumed 
the measurable broadening of profiles seen in the C I and K I spectra, beyond the instrumental smoothing, 
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was created by clusters of narrow profiles rather than a smooth turbulent broadening. This conservative 
model indicates that as the value of /A increases by one order of magnitude from 5 to 50, which corresponds 
to peak values of t(v) equal to 0.24 and 2.4, respectively, the average discrepancy in r a //A over the strongest 
portions of the profile (—36 < v < —28 km s" 1 ) is only 0.041 dex. This is only about 13% of the observed 
rate of divergence between our /-values and those that have been published elsewhere. 

We have demonstrated the expected insignificance of misrepresentations reflected by T a when progressing 
from very weak absorptions to absorptions of moderate strength (10 times stronger), where we consider r a to 
be simply an instrumentally smoothed version the real optical depth r. Moving on, we are now in a position 
to interpret comparisons between the observed absorption features and ones that are reconstructed from our 
derived column densities vs. velocity, on the assumption that the computed values of r a give the correct 
residual intensities (to within about 0.041 dex). The different displays in Fig. [16] show this comparison 
for the complex at —34 km s _1 for both the /-values derived by JT01 and those listed by Froese Fischer 
& Tachiev (2004) (FT04). We have restricted the choices in this display to transitions of C I*, and we 
have rejected cases where there was interference from other transitions that had absorptions from other, 
less negative velocity components that were overlapping. It is clear that as the lines become weaker, the 
disparity between the observations and the predictions using the FT04 /-values increases, while those using 
the JT01 /-values match the observations to within the noise levels in each case. 

It is important to emphasize that the investigation discussed here addresses the accuracy of the relative 
/-values from one multiplet to the next and not the correctness of all of the multiplets taken together. Had 
we adopted a multiplet other than the strongest one at 1657 A as a standard, all of our derived /-values 
would have been lower by some constant factor. 
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Fig. 16. — Selected absorptions from the —34 km s _1 C I* component in the spectrum of HD 210839 (A Cep) 
that arc not contaminated by overlapping features created by the components between —20 and Okm s _1 . 
The observed fluxes (histogram- style traces) have been normalized to our best estimate for the continuum 
level in each case. The different panels are ordered in a sequence of decreasing values of log(/A) determined 
by JT01. The upper smooth curve in each panel represents a reconstruction of the profile using the published 
f-value, while the lower curve arises from the f-value derived by JT01. [Key to published /-value sources: 
M03 = Morton (2003), FT04 = Froese Fischer & Tachiev (2004).] 
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